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EXECUTIVE  SUMMARY 


Current  vulnerability/lethality  analysis  methodology  is  concerned  strictly  with 
the  local  damage  caused  by  penetration  of  KE  penetrators  or  shaped  charge  jets. 

These  damage  mechanisms  arc  of  primary  interest  and  have  been  extensively  studied, 
both  from  a  theoretical  and  a  practical  basis;  however,  another  damage  mechanism  of 
considerable  significance  is  ballistic  shock.  As  tactical  combat  vehicles  are  becoming 
better  protected  against  perforating  impacts,  the  incidence  of  component  and  structure 
failure  due  to  ballistic  shock  is  becoming  more  significant  in  terms  of  the  overall 
damage  levels.  Current  methodology  does  not  predict  or  account  for  these  eff  :'ts  in  a 
satisfactory  and  complete  manner.  This  report  will  attempt  to  address  one  a  ct  of 
this  complex  problem  as  a  genesis  to  the  development  of  a  suitable  vulnerability/ 
lethality  methodology  that  includes  ballistic  shock  damage  mechanisms. 

Whereas  the  damage  directly  associated  with  the  penetration  event  is  considered 
local  in  nature,  there  can  be  global  effects  that  degrade  and  incapacitate  vehicle 
systems.  These  global  effects  arise  from  the  forced  vibratory  loading  due  to  the 
impact  of  the  attacking  munition.  Such  effects  include  far  field  structural  failures, 
inoperable  electro-optical  devices  and  biological  damage.  As  an  initial  assessment, 
this  report  will  consider  the  gross  motion  of  heavy  structures  subjected  to 
nonperforating  impacts  (impulse  loads)  and  what  effect  the  structure  interface  has  on 
the  system  dynamics.  The  structures  are  modeled  as  a  single-degree-of-freedom 
(SDOF)  system  and  a  two-degree-of-freedom  (2DOF)  system  to  represent  a  turretless 
vehicle  and  a  turreted  vehicle,  respectively.  The  turret  and  hull  are  idealized  as  rigid 
bodies  connected  together  by  a  spring  interface  of  variable  stiffness.  The  dynamical 
equations  of  motion  arc  derived  under  the  condition  of  no  slip  and  consist  of  two 
independent  equations  for  each  of  the  four  main  regimes  of  motion  that  the  2DOF 
system  can  occupy.  The  second  order,  ordinary  differential  equations  are  coupled  and 
nonlinear.  Analytic  solutions  are  determined  by  linearizing  and  uncoupling  the 
equations  via  modal  analysis  techniques.  Due  to  physical  discontinuities  at  the  impact 
points,  these  analytical  solutions  are  only  valid  in  the  continuous  regimes  of  motion  of 
the  2DOF  system.  Conservation  of  moment  of  momentum  is  used  to  derive  transition 
equations  that  bridge  these  discontinuities  by  allowing  the  calculation  of  new  initial 
conditions  for  the  next  regime  of  motion  that  the  system  enters.  The  local  time 
histories  of  each  regime  of  motion  are  then  pieced  together  to  yield  the  global 
dynamical  response.  By  incorporating  these  analytical  solutions  within  a  FORTRAN 
computer  simulation,  the  gross  motion  of  the  vehicles  under  varying  impacts  can  be 
determined.  A  parametric  study  of  the  motion  is  performed  by  varying  the  stiffness  of 
the  turret  interface  and  the  physical  properties  of  both  the  hull  and  turret.  It  is  shown 
that  system  stability  is  a  strong  function  of  the  mass  and  geometry  and  a  weak  function 
of  the  stiffness  parameter.  It  is  also  shown  that  the  mass  of  the  turret  has  a  strong 
damping  effect  on  the  motion  of  the  hull  due  to  the  spring  interface.  This  effect  is 
more  pronounced  the  stiffer  the  interface  becomes.  In  general,  it  can  be  concluded  that 
though  the  system  stability  is  highly  dependent  on  the  mass  and  geometry  of  both  the 
turret  and  hull,  the  effect  of  the  stiffness  associated  with  the  turret  interface  can 
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stabilize  the  system  under  ceratin  conditions  and  as  the  stiffness  of  the  turret  interface 
increases,  the  response  of  the  2DOF  system  becomes  more  like  a  rigid  body,  which  is 
to  be  expected. 
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1.  INTRODUCTION 


With  the  advent  of  live-fire  testing  involving  heavily  armored  structures,  it  has 
become  apparent  that  damage  mechanisms  other  than  penetration  and  perforation  have 
considerable  significance  in  terms  of  munition  lethality  and  taiget  vulnerability. 
Specifically,  ballistic  shock  effects  have  become  more  pronounced  in  terms  of  the 
overall  damage  levels  sustained  by  armored  targets  during  these  tests,  both  for 
perforating  and  nonperforating  impacts.  Though  the  taiget  damage  assessment 
includes  the  effects  associated  with  ballistic  shock,  current  vulnerability  and  lethality 
analyses  do  not.  In  light  of  this  situation,  the  U.S.  Army  has  identified  the  need  to 
better  understand  the  ballistic  shock  phenomenon  and  to  modify  and  upgrade  current 
vunerability  and  lethality  methodologies  to  satisfactorily  address  this  issue. 

An  in-house  research  effort  was  initiated  during  the  1991  time  frame  by  the  U.S. 
Army  Ballistic  Research  Laboratoiy  (BRL),  which  was  transitioned  into  the  current 
U.S.  Army  Research  Laboratory  (ARL).  The  programmatic  approach  was  to  look  at 
various  ways  to  calculate  and  predict  damage  levels  associated  with  ballistic  shock 
and  incorporate  the  most  reasonable  method  within  a  current  vulnerability 
methodology,  with  reasonable  being  defined  as  resource  efficient,  technically  correct 
and  answers  that  arc  simple  and  robust,  so  to  speak.  This  report  details  the  research 
efforts  of  one  of  the  methods  proposed  to  address  certain  areas  within  the  overall 
ballistic  shock  program,  the  rigid  body  method  (RBM).  The  author  believes  that  this 
method  has  its  greatest  utility  in  addressing  ballistic  shock  effects  in  terms  of  crew 
casualty  predictions  and  assessments  due  to  acceleration  levels  and  structural  failure 
of  the  turret-hull  interface.  Before  discussing  the  research  that  forms  the  basis  of  the 
method  in  detail,  it  is  appropriate  to  discuss  shock  in  general  and  the  applicability  of 
this  method. 

Shock  is  generally  defined  as  a  relatively  laige  force  applied  suddenly  and 
quickly  with  a  time  period  that  is  relatively  short  as  compared  to  the  natural  period  of 
the  structure  that  is  being  subjected  to  this  force.  This  transient  force  can  produce 
damage  local  to  the  point  of  application  of  this  force  and  also  vibratory  forces  that 
affect  the  structure  beyond  the  local  point  of  application.  These  nonlocal  effects  are 
called  the  “global  effects”  and  they  are  the  focal  point  of  this  research  effort. 

The  analytical  technique  termed  “RBM”  was  investigated  for  a  number  of 
reasons.  First  and  foremost,  the  assumptions  associated  with  this  method  allow  the  use 
of  some  powerful  mathematical  techniques  such  as  modal  analysis,  that,  in 
conjunction  with  engineering  dynamics,  provide  well-formulated  equations  of  motion 
that  are  solvable  either  analytically  or  numerically.  Secondly,  this  method  should 
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provide  the  most  conservative  answer  of  the  proposed  methods  in  terms  of  structural 
survivability,  thus  setting  an  upper  bound.  Lastly,  modeling  of  the  physical  system 
becomes  relatively  simple  and  unambiguous,  requiring  only  computational  algorithms 
and  methods  that  are  not  time  intensive.  For  these  reasons,  this  method  was 
investigated  and  the  results  generated  using  this  method  will  be  compared  with 
empirical  data  generated  from  test  programs  when  available. 

A  brief  technical  synopsis  of  this  method  and  general  vibration  theory  is 
presented  here  to  facilitate  a  greater  understanding  and  to  present  the  advantages  and 
limitations  of  this  approach,  as  seen  at  this  time.  A  structure,  whether  considered  as  a 
single  component  with  isotropic  material  properties  or  as  a  conglomeration  of 
components  that  may  or  may  not  share  similar  material  properties,  has  associated  with 
it  a  structural  parameter  called  the  natural  frequency(s).  The  natural  frequency,  being  a 
function  of  the  material,  geometry,  and  boundary  conditions,  is  an  inherent  property  of 
the  structural  system  only,  not  any  external  conditions.  Thus  the  natural  frequency,  or 
free  vibration  frequency,  is  the  frequency  or  frequencies  at  which  the  structure 
oscillates  after  the  external  forcing  function,  i.e,,  the  initiating  impulse,  is  removed. 
This  natural  frequency  (or  frequencies)  of  a  structure  define  the  structure’s  response  to 
external  forces;  i.e.,  when  the  structure  is  subjected  to  a  broad-frequency  acceleration 
environment  such  as  an  impact,  the  structure  will  absorb  the  energy  more  easily  at 
certain  frequencies,  the  natural  frequencies  of  the  system.  So  the  acceleration  response 
spectrum  of  the  shock  environment  will  have  specific  values  at  each  natural  frequency, 
the  sum  of  which  determines  the  total  response  of  the  structure.  These  acceleration 
values  at  each  natural  frequency  specify  the  deformations  and  stresses  induced  in  the 
structure  with  the  lowest  frequency  referred  to  as  the  fundamental  frequency.  The 
greatest  deformation,  and  therefore,  the  largest  stresses  and  strains,  occur  at  the 
fundamental  frequency  of  the  structure,  in  general.  It  is  this  fundamental  frequency 
that  we  are  interested  in  because  it  is  associated  with  the  largest  peak  acceleration 
values,  and  therefore,  has  the  greatest  potential  for  damage.  If  we  can  assume  that  the 
fundamental  frequency  of  a  real  structure  is  close  in  magnitude  to  a  similar  structure 
modeled  as  a  rigid-body,  then  this  method  should  be  useful. 

This  idealized  assumption  that  the  structure  behaves  as  a  rigid  body  requires 
that  the  structural  stiffness  or  rigidity  be  high,  or  conversely,  the  structural  damping 
approach  zero  since  the  damped  frequency  is  defined  as  the  natural  frequency  of  a 
system  with  damping,  or  mathematically: 

wd.wn«JiT^  (i) 

where  £  -  c  /  (2  m  Wn ) 

Wd  -  damped  frequency 

Wn  -  natural  frequency  of  system  with  no  damping 
m  -  mass 

c  -  damping  constant 
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Thus,  when  the  damping  constant  c  approaches  zero,  then  the  damped  frequency 
approaches  the  natural  frequency,  as  shown  in  equation  1.  If  we  assume  a  structure 
composed  of  a  very  small  number  of  components,  we  could  determine  the  natural 
frequencies  of  each  of  these  simple  components,  but  we  could  not,  in  general, 
determine  the  natural  frequencies  of  the  structure  from  these  individual  frequencies. 
The  structural  frequencies  are  functions  of  the  mass,  stiffness,  and  damping  of  each 
component  and  the  stiffness  and  damping  asssociated  with  the  interface  between 
components.  Thus  to  analyze  anything  but  an  extremely  simple  model  requires  a  finite 
element  analysis  (FEA)  program.  So  the  following  assumptions  are  necessary  for 
rigid-body  motion.  The  ratio  of  the  stiffness  of  each  individual  component  and  their 
interfaces  to  the  stiffness  of  an  equivalent  structure  of  homogeneous  material  should 
approach  unity,  where  the  stiffness  is  a  function  of  material,  material  impedance  at  the 
interfaces,  and  structural  rigidity.  Thus  the  damping  terms  must  approach  zero  or  be 
very  small  compared  to  the  total  mass,  as  shown  in  equation  1.  For  an  application  of 
very  heavy  targets  with  rigidly  welded  plates,  this  assumption  is  considered  reasonable 
at  this  time. 

Assuming  rigid-body  motion,  we  can  model  an  armored  vehicle  as  a  simple 
arrangement  of  two  blocks  connected  by  an  equivalent  spring  constant  representing 
the  stiffness  of  the  bolted  turret/hull  Interface.  This  is  a  two-degree-of-freedom 
(2DOF)  system  requiring  the  solution  of  two  simultaneous,  second  order,  nonlinear, 
inertially  and  elastically  coupled  equations.  This  system  can  be  solved  two  ways, 
either  numerically  or  analytically.  Analytic  solutions  for  this  problem  were 
determined;  though,  in  general  there  are  very  few  known  analytic  solutions  to 
nonlinear  equations.  In  order  to  do  this  the  nonlinear  equations  were  linearized  and 
then  uncoupled  using  modal  analysis  techniques.  This  procedure  allows  for  the 
relatively  simple  generation  of  analytic  solutions  for  various  forcing  functions  where 
the  forcing  function  chosen  is  an  analytic  representation  or  model  of  the  shock 
producing  impact.  The  analytic  representations  under  consideration  for  the  forcing 
function  arc  the  Dirac  delta  function  (actually  a  distribution)  and  a  single  pulse  sine 
wave.  Both  of  these  functions  were  examined  for  derivation  of  the  equations  of  motion 
but  the  system  examined  was  the  free  vibration  case  since  the  purpose  of  this  initial 
report  is  methodology  development  only.  Development  and  validation  of  the  specific 
forcing  functions  for  ballistic  shock  modeling  is  currently  being  accomplished  by  the 
University  of  Dayton  Research  Institute  (UDRI)  as  part  of  the  overall  shock  program. 
Finally,  one  of  the  long  term  goals  would  be  to  solve  the  system  of  nonlinear  equations 
numerically  for  comparison  to  the  linear  solutions. 

The  majority  of  research  efforts  dealing  with  rigid  body  motion  lies  in  the  field 
of  civil/mechanical  engineering,  specifically,  structural  dynamics.  A  literature  search 
was  conducted  and,  though  the  relevant  papers  had  strong  civil  orientations,  the 
underlying  physical  and  mathematical  principles  apply  to  this  particularly  military 
application. 

The  prevailing  trend  in  structural  dynamics  is  to  assume  that  structures  subjected 
to  dynamic  loads  or  excitations  behave  as  a  deformable  continuum.  There  are 
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numerous  examples  of  physical  phenomena  that,  under  certain  conditions,  contradict 
this  assumption  and,  in  fact,  show  that  rigid  body  motion  is  a  significant  response 
mechanism.  The  collapse  of  the  Cypress  Street  Viaduct  (1880)  during  the  1989  Loma 
Pricta  earthquake  was  kinematic  in  nature  [1]  and  the  Transamerica  building  exhibited 
a  significant  amount  of  rotation  or  rocking  about  its  base  [2]  during  the  same  event. 
Other  examples  of  this  rigid  body  behavior  include  the  cyclic  rocking,  with  the 
corresponding  stretching  of  the  base  anchor  bolts,  of  the  relatively  tall  and  slender 
petroleum  cracking  towers  during  the  Arvin-lbhachapi  earthquake,  California,  21  July 
1952  [3],  The  primary  reason  for  loss  of  structural  integrity  in  these  events  was  the 
rocking,  and  in  some  cases,  subsequent  toppling  of  these  structures,  and  not  due 
explicitly  to  some  material  failure  mechanism.  These  events  highlight  the  need  for  a 
better  understanding  of  how  structures  undergoing  rigid  body  motion  respond  to 
external  stimuli. 

The  simplest  way  to  model  this  behavior  is  by  analyzing  the  rocking  motion  of  a 
single  rigid  block  on  a  rigid  foundation.  This  has  been  done  extensively  in  the 
literature,  in  recent  times  by  Housner  [3]  and  then  expanded  upon  with  variations  to 
include  flexible  foundations  of  the  two  spring  type  [4]  and  the  Winkler  type  [4,5,6], 
and  both  harmonic  and  random  base  excitations  [7,8,9,10].  A  comprehensive  review 
of  the  research  up  to  1980  was  compiled  by  Ishiyama  [11]  and  is  a  good  introduction 
to  this  field.  A  more  recent  paper  by  Lipscombe  and  Pellegrino  [12]  compares 
experimental  data  from  single  rocking  blocks  to  predicted  data  from  the  theoretical 
models  of  such  systems  with  some  interesting  insights  and  conclusions. 

This  characterization  of  a  structure  as  a  single  rigid  body,  though  desirable  in 
terms  of  its  simplicity,  does  not  adequately  address  the  complexities  inherent  in  most 
typical  man-made  structures  of  interest.  Thus  a  few  researchers  have  extended  the 
research  from  single-degree-of-frcedom  (SDOF)  systems  such  as  the  single  rocking 
block  to  multi-degree-of-frcedom  (MDOF)  mechanisms,  some  as  simple  as  two 
stacked  blocks  [13]  and  others  as  complex  as  multiple  story  structures  composed  of 
individual  rigid  bodies  [14,15,16,17].  These  researchers  have  considered  systems  of 
interest  comprised  of  rigid  bodies  stacked  upon  each  other  with  no  provision  for  any 
type  of  fastening  arrangement.  These  rigid  body  assemblies  represent  precast  concrete 
building  systems  without  considering  any  physical  connection  between  the  individual 
rigid  bodies.  This  report  will  attempt  to  extend  the  understanding  of  the  behavior  of 
such  structures  by  addressing  the  dynamics  of  a  rigid  body  system  whereby  the  rigid 
bodies  are  interconnected  or  fastened  together  with  springs  of  variable  resistance  to 
model  physical  connectors  such  as  reinforcing  bars,  bolted  connections,  or  a  turret/ 
hull  interface.  This  report  begins  to  address  what  effect  the  fasteners,  geometry 
variations,  and  mass  considerations  have  on  the  system  dynamics. 

« 

One  idealization  of  such  a  MDOF  system  can  be  made  by  considering  two 
symmetrically  stacked  blocks  interconnected  by  two  springs  where  the  top  block  can 
represent  the  turret  of  an  armored  vehicle  with  the  bottom  block  representing  the  hull. 
The  blocks  arc  allowed  to  rotate  about  their  edges  and  the  coefficient  of  friction  is 
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considered  large  enough  to  preclude  any  slip;  thus,  the  poles  of  rotation  can  be  thought 
of,  in  a  kinematic  sense,  as  a  pin  constraint.  A  mathematical  model  is  developed  for 
such  a  system  and  the  resulting  equations  of  motion  are  solved  analytically.  The 
solutions  are  encoded  within  a  computer  program  and  the  time  histories  of  the  two 
blocks  are  examined  with  respect  to  parameter  variations  of  the  masses  and  spring 
constants. 

Since  the  time  histories  are  determined  as  functions  of  the  absolute  rotation 
angles  of  each  block,  and  kinematic  expressions  can  be  derived  for  velocities  and 
accelerations  at  any  point  of  the  rigid  bodies  as  a  function  of  these  absolute 
coordinates,  then  the  g-lcvels  at  any  point  can  be  determined  to  assess  crew 
incapacitation  values.  Also,  failure  of  the  bolted  interface  between  the  hull  and  turret 
can  be  calculated  for  known  material  and  geometry  specifications.  The  modeling  and 
derivations  are  presented  next. 
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2.  MATHEMATICAL  MODELING 

The  first  system  to  consider  for  simplicity  and  possible  insights  into  more 
complex  systems,  is  a  single  rigid  block  rocking  on  a  horizontal  surface  with  attached 
springs.  The  assumptions  with  this  model  are  that  the  springs  operate  only  in  the 
Linear  elastic  range  and  are  massless.  See  Figure  1  below. 


Y 

k 


i 

i 

. 

2, 

\ 

1 

1 

I 

Pi  J  t 

'•  ••  'V.'  ::::i 

1 

P» 

/ /  A 

H 

_i _ ! 

7  77 
R; 

r 

y  /  A 

0  \ 

7  7777 

r 

2j 

b; - ► 
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2  bj  “->  Width  of  Block  Bj  hj  — >  Height  of  Center  of  Mass 

C  ]  Center  of  Mass  n  j  >  Distance  from  Edge  to  Center  of  Mass 

l],  I2  ■•>  Distance  from  Block  Center  to  Springs  Ri  and  R2  respectively 

FIGURE  1  -  Single  Block  with  Springs  R]  and  R2 


Now  consider  block  B\  with  some  initial  angular  displacement,  as  shown  in 
Figure  2.  Note  that  the  block  rocks  about  its  comers  P\  and  P2,  and,  for  the  shown 
configuration,  spring  R[  remains  at  a  constant  length  while  spring  R2  is  in  tension. 
The  block  will  eventually  rotate  clockwise  and  impact  the  horizontal  surface 
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Cj  -—>  Center  of  Mass 
H  —  *  >  Length  of  Springs  Rj,  R2 

K},K2  — - >  iSpr/n#  Constants  of  Springs  Rj  and  R2  respectively 
L  - — >  Length  of  Springs  Rj  andR2>  as  an  implicit  function  ofQ 

FIGURE  2  -  Single  Block  with  Displacement  0 

and  then  begin  to  rotate  about  point  P2.  It  is  this  transition,  from  rocking  about 
comer  P\  to  rocking  about  comer  P2,  that  is  one  source  of  nonlinearity  for  this 
seemingly  simple  problem.  Thus  the  equation  of  motion  describing  this  rocking  will 
be  discontinuous  at  8  =  0,  which  implies  that  the  character  of  the  equation  will 
differ  according  to  the  sign  of  the  angle  of  displacement,  6.  Taking  moments  about 
points  P[  and  P2  gives  the  following  equations  of  motion. 


Ipi  0  +  Mg  (n jcosb-h }Sin9)+K2(  1  -HH)( b j +l2)[Hcosb+(b j+l2 )sinQ\-0  (2a) 

for  0  >  0  where  Ipj  « I  a  +  M(nj2  +  h2) 


Ip2$  •  Mg[(2  b]-nj)cosQ+hjSinQ))+Kj(l-H/L)(b]+lj)[HcosB+(hj+lji)sinQ]=0  (2b) 
for  0<O  where  Ip2  =  I&  +  M[(2 bj-nj?+ hf] 


These  are  the  general  forms  of  the  equations  of  motion.  If  one  assumes  small 
displacements  (  0  <  20°  ),  a  homogeneous  distribution  of  mass,  and  that  no  gap 
exists  between  the  ground  and  the  reference  frame  Ctf-0),  then  the  simplified 
equations  are: 
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hi  6  +  [Krfbi+hf-Mgaj]  0  =  -Mgbj 
In  ®  +  [Krfbj+ljf-Mgaj]  9  =  A#**; 


for 

0>O 

(3a) 

for 

e<o 

(3b) 

An  explanation  concerning  the  validity  of  the  previous  assumptions  is  required  at 
this  time.  One  assumes  small  displacements  in  order  to  eliminate  the  nonlinearities 
associated  with  the  trigonometric  functions.  Thus,  the  simplified  equations  are  valid 
only  as  long  as  the  rotation  of  the  rocking  block  remains  within  the  range  of  the  small 
angular  displacements.  The  second  assumption  concerning  a  homogeneous 
distribution  of  mass  simply  implies  that  the  block’s  center  of  gravity  coincides  with 
the  geometric  center.  This  is  a  reasonable  assumption  for  most  man-made  structures 
that  would  be  considered  as  rigid  body  structures  exhibiting  the  type  of  response  this 
study  addresses.  The  last  assumption  is  more  subtle  than  the  other  two  and  requires  a 
more  detailed  explanation. 

As  stated  earlier  in  this  report,  H  represents  the  free  length  of  the  spring,  but  in 
actuality,  it  is  only  an  artifact  of  how  the  system  has  been  modeled .  H  could  have  been 
modeled  as  the  distance  from  the  horizontal  base  to  the  block,  i.e.,  the  increase  in  the 
length  of  the  spring,  with  a  spring  free  length  of  zero.  This  would  effectively  remove 
H  from  the  equations  since  the  spring  force  acting  on  the  block  is  a  function  of  6  only. 

There  are  two  important  aspects  or  features  that  equations  3a  and  3b  bring  to  light 
for  this  relatively  simple  system  which  should  hold  true  for  more  complex  or  MDOF 
systems  also.  The  right-hand  side  (RHS)  changes  sign  from  equation  3a  to  equation 
3b.  This  nonlinearity,  mentioned  previously,  is  due  to  the  transition  front  rocking 
about  one  comer  to  the  other  comer.  The  second  concerns  the  terms  or  coefficients 
associated  with  the  displacement  term,  0.  The  addition  of  the  springs  to  the  problem 
can  change  the  character  of  the  equations.  If  the  spring  constants  K\  and  K2  we  set 
equal  to  zero,  then  the  equations  have  negative  stiffness.  If  the  spring  constants  and 
certain  physical  parameters  are  great  enough  in  magnitude,  then  it  is  possible  for  the 
equations  to  have  positive  stiffness.  This  is  a  significant  feature  which  directly  affects 
the  form  of  the  solution  to  these  equations,  as  will  be  discussed  next. 

The  solution  set  to  equations  3a  and  3b  is  presented  below  in  equations  4a  and 
4b,  respectively.  Since  the  equations  of  motion  are  discontinuous  at  0  -  0,  a  solution 
set  is  provided  for  both  conditions,  0  >  0  and  0  <  0  by  using  the  sgn  function. 


0(f)  =  Dj  cos  JA  t  +  D2  sin  JA  t-  (B/A)  sgn  0 
0(f) «  Dj  cosh  jAt  +  D2  sinh  JA  t  +  (B/A)  sgn  0 


for 

A>  0 

(4a) 

for 

A<  0 

(4b) 
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where  A  =  ]A|  since  the  sign  of  A  is  accounted  for  hy  the  form  of  the  solutions 
and  specifically:  A  =  [K2(bj+I2f  -  MgajJ  /  Ipi  =  [ K}(bj+l])t2  -  Mgaj]  1  b2 
B  =  Mgbflpj  ~  Mgbj/lp2 

Dj-  Q(0)  +  BIA  D2  =  Q(0)iJA  for  A>0 

Dr-*tyO)-BtA  D2  =  b(0)/jA  for  A  <0 


Thus  it  is  evident  that  the  sign  of  the  displacement  angle  0  affects  only  the  particular 
term  of  the  total  solution  while  the  sign  of  the  parameter  A  affects  both  the 
homogeneous  part  and  the  particular  part  of  the  total  solution.  The  significance  of  this 
is  that  once  the  geometry  of  the  block  is  defined  then  the  sign  of  parameter  A  is 
determined  by  the  value  of  the  stiffness  term  K /.  The  sign  of  A  is  determined  by 
equation  Sa  and  5b  presented  below  and  as  initially  defined  above. 


A  «  [K2(bi+l2?  -  Mgaj]  /  lP1 
A  =  [K1(b1+l1jl-MgaII/Ip2 


for  0  >  0 

(5a) 

for  0  <  0 

(5b) 

By  setting  the  parameter  A  to  zero  in  equations  5a  and  5b  and  solving  for  the  Kt  term, 
the  threshold  value  for  when  the  character  of  the  solution  changes  from  equation  4a  to 
4b  or  vice  versa  can  be  determined  readily.  These  equations  are  presented  next. 


K2  =  Mgaftbj+li? 

for  e>o 

(6a) 

K,  =  MgajHbj+lj? 

for  8  <0 

(6b) 

Titus,  when  the  value  of  K2  is  greater  than  the  value  as  determined  in  equation  6a  then 
the  solution  is  sinusoidal  in  nature,  and,  if  less,  the  solution  is  hyperbolic  in  character, 
A  similar  statement  can  be  made  for  the  Kj  term.  The  practical  significance  of  this 
phenomenon  will  be  discussed  in  the  results  section  of  this  report. 

The  nonlinearity  due  to  transition  of  the  block  from  the  0  >  0  case  to  the  0  <  0 
case  requires  that  the  solution  consist  of  two  solutions,  one  from  each  of  the  linear 
regimes.  Thus  the  total  or  global  solution  will  consist  of  two  local  solutions  which  can 
be  written  in  closed  form  due  to  the  linear  nature  of  the  equations.  This  leads  to  an 
important  consideration:  how  does  the  transition  from  rocking  about  one  comer  to 
rocking  about  the  other  comer  affect  the  global  solution?  The  two  local  solutions  are 
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tied  together  by  simply  using  the  final  conditions  of  the  pre-impact  or  pre-transition 
regime  as  the  initial  conditions  for  the  next  regime.  This  implies  that  there  is  no 
energy  loss  during  the  transition  and  one  would  expect  that  the  block  would  rock 
through  an  angle  9  -  -0O  where  0O  the  initial  displacement  of  the  block  (assuming  no 
initial  velocity).  Therefore,  the  block  would  continue  to  rock  forever  for  this 
conservative  case  of  no  energy  dissipation.  Practically,  this  case  is  not  tr  ue,  because 
during  the  transition,  which  is  really  an  impact  of  the  block  with  its  supporting 
surface,  there  is  an  energy  loss.  In  classical  mechanics,  this  energy  loss  is  accounted 
for  by  a  proportionality  constant  (more  formally  known  as  a  restitution  coefficient,  C,) 
that  relates  the  pre-impact  velocity  to  the  post-impact  velocity  such  as  equation  7. 


6f-C$i  where  0£C<1  (7) 


Equating  the  kinetic  energy  of  the  block  immediately  prior  to  impact  and  immediately 
after  impact  (since  the  potential  energies  are  equivalent)  with  a 
proportionality  factor,  /?,  to  account  for  the  reduction  in  energy  due  to  the  impact 
gives  equation  8. 


(8) 


Realizing  that  IP1  =  IP2  for  a  homogeneous  block,  then  the  reduction  in  energy  due  to 
impact  can  be  represented  as  R  below. 


R- 


$ 


6 


t 

2 


i.e.,  R  *  C2 


(9) 


Now,  if  the  impact  can  be  considered  as  perfectly  plastic,  i.e.,  the  block  does  not 
bounce  or  slip,  then  the  block’s  rotation  is  smooth  and  transits  to  the  other  pole  of 
rotation  in  such  a  manner  that  conservation  of  angular  momentum  can  be  assumed. 
This  assumption  is  pivotal  in  that  it  allows  the  derivation  of  an  analytical  expression 
relating  the  pre-impact  and  post-impact  velocities  to  account  for  the  discontinuity  that 
occurs  at  0  =  0.  See  the  appendix  for  the  derivation  of  the  general  equation  for 
conservation  of  moment  of  momentum.  Equating  the  moment  of  momentum  about 
point  P[  immediately  prior  to  impact  and  immediately  after  impact  results  in  the 
following  relationship. 


IPj  $1  -2  mb]  &i(b}  cosQ  -  aj  sinQ)  -  IP2&f  noting  lPj  ~  lP2  (10) 
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And  since  0  **  0  during  impact,  equation  10  can  be  simplified  with  some  rearranging. 


C- 
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\a]  +  b]  > 
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What  this  equation  implies  is  that  even  if  one  assumes  the  idealized  conditions  of  a 
smooth  transition  from  one  pole  of  rotation  to  the  other  and  rigid  body  conditions, 
there  is  going  to  be  an  energy  loss  during  impact  that  is  strictly  a  function  of  the 
geometry  of  the  block!  Thus,  using  equation  9,  this  energy  loss  can  be  represented  as 


energy  loss  -  /  -  R  or  energy  loss  -  1  -  C2 


(12) 


To  minimize  the  energy  loss  during  transition  due  to  the  block  geometry,  the  last  term 
of  equation  11  must  approach  zero,  which  mathematically  implies  that  a;  >>  bj,  i.e., 
tall,  slender  blocks.  Conversely,  short,  wide  blocks  where  bj»  aj  corresponds  to  the 
maximum  amount  of  energy  loss.  Recognizing  that  the  theoretical  lower  limit  of  C  is 
zero,  it  is  a  trivial  exercise  to  show  that  this  occurs  when  the  ratio  of  b\  to  a\ 
approaches  the  value  1.414  or  Ji.  Also,  it  is  interesting  to  note  that  if  the  ratio  of  b}  to 
a i  is  less  than  1.414,  then  C  will  be  negative!  Mathematically,  this  implies  that  the 
post-impact  velocity  is  of  opposite  sign  from  the  pre-impact  velocity  so  there  is  a 
instantaneous  change  it;  the  direction  of  rotation  of  the  block.  Physically,  this  means 
that  the  pole  of  rotation  of  the  block  does  not  change  to  the  other  comer,  but  that  the 
block  bounces  back,  maintaining  the  same  pole  of  rotation.  Practically,  this  does  not 
occur  and  this  condition  implies  that  the  block  bounces,  slips  and  leaves  the  surface  of 
contact,  thereby  changing  the  scope  of  this  problem  and  invalidating  the  assumption  of 
conservation  of  moment  of  momentum.  A  recent  paper  by  Llpscombc  and  Pellegrino 
[12]  discusses  this  phenomenon  in  greater  detail  and  concludes  that  the  SDOF  model 
with  the  assumption  of  conservation  of  moment  of  momentum  is  valid  for  slender 
blocks,  but  that  this  method  is  often  inaccurate  for  short  blocks  where  bouncing  tends 
to  be  more  pronounced.  The  paper  also  presents  a  comparison  of  experimental  data 
for  varying  block  geometries  and  three  methods  to  model  the  phenomenon  termed  as 
“bouncing”. 

Although  the  condition  represented  by  equation  6  with  the  theoretical  limits  of 
the  proportionality  factor  C  to  model  the  impact  or  transition  is  necessary,  it  is  not 
sufficient.  By  using  the  conservation  of  angular  momentum  principle,  a  more 
conservative  estimate  of  the  theoretical  limits  of  C  based  on  the  geometry  of  the  rigid 
body  can  be  determined.  This  principle  will  also  apply  to  more  complex  cases,  as  is 
discussed  next. 

The  next  step  is  to  consider  a  more  complex  system  or  MDOF  system  (Figure  3), 
now  that  the  relatively  simpler  system  has  provided  some  insights  into  the  modeling 
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2  aj  -->  Height  of  Block  Bj  lj,  l2  -->  Distance  from  Block  Center  to  Spring 

?  bj  ~>  Width  of  Block  Bj  hph2  ->Height  of  Center  of  Mass 

2  a2  - >Height  of  Block  B2  n  j  n2-> Distance  from  Edge  to  Center  of  Mass 

2  b2  ->  Width  of  Block  B2  r  j  ft  2-->  Springs  Rj  and  R2 

C}~>Mass  Center  Block  B /  C2-->Mass  Center  Block  B2 

FIGURE  3  -  Model  of  Two  Stacked  Blocks 

and  form  of  the  equations  of  motion.  This  system  consists  of  two  blocks  stacked 
together  that  are  connected  via  springs  and  arc  placed  on  a  rigid  horizontal  surface. 
All  contact  conditions  are  no-slip,  the  springs  are  massless  and  linear,  and  the  system 
is  subject  to  free  vibration  only,  i.e„  no  external  applied  forces.  The  model  of  the 
system  portrays  the  two  blocks  as  separated  by  the  length  of  the  two  springs  but  in 
actuality  the  system  under  consideration  in  this  report  has  the  biocks  in  direct  contact. 
This  was  just  done  as  a  visualization  aid  during  the  modeling  process.  As  with  the 
single  block  model  where  this  distance  was  represented  by  H,  the  same  modeling 
artifice  was  used  here  and  it  is  assumed  that  H  is  zero  to  get  the  cbrrect  equations  for 
this  particular  model. 

This  system  of  two  rocking  blocks  exhibits  the  same  character  as  the  single 
rocking  block  in  the  sense  that  the  equations  change  when  the  displacement  angles 
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change  sign.  This  follows  from  the  mathematical  and  physical  discontinuity  that 
occurs  whenever  either  block  transitions  from  rocking  about  one  comer  to  rocking 
about  the  other  comer.  In  addition,  the  point  about  which  the  top  block  rocks  will  not, 
i'  general,  be  stationary  or  fixed  as  it  is  for  the  bottom  block;  thus,  there  are  Coriolis 
terms  introduced  into  the  final  equations  of  motion.  Furthermore,  the  motion  of  each 
block  is  highly  dependent  or  coupled  to  die  motion  of  the  other  block  introducing 
further  complexity  into  a  seemingly  simple  problem.  The  aggregate  effect  of  these 
system  characteristics  results  in  equations  of  motion  that  are  highly  nonlinear  and 
atypical. 

To  make  this  problem  more  tractable,  at  least  from  a  mathematical  viewpoint,  one 
must  differentiate  between  the  various  regimes  of  motion  that  this  system  can  possibly 
attain.  Each  regime  of  motion  can  be  defined  by  the  particular  equation  of  motion  that 
distinguishes  one  regime  from  another.  It  is  somewhat  obvious,  especially  considering 
the  insights  gained  by  the  single  rocking  block  problem,  that  this  system  can  rock  in  at 
least  four  modes  or  phases,  based  on  the  geometry  of  the  system  and  how  the 
displacement  angle  sign  changes  occur.  There  are  two  other  subphases  that  will  be 
considered  only  as  special  cases  of  these  four  general  phases.  The  four  phases  are 
presented  in  Figure  4. 

The  two  special  cases  occur  when  both  blocks  rock  in  tandem  (rigid  body  mode) 
and  when  only  the  top  block  rocks,  which  is  simply  the  problem  of  the  single  rocking 
block  discussed  earlier.  The  only  interest  in  the  special  case  of  both  blocks  rocking  in 
concert,  which  once  again  can  be  considered  as  the  single  rocking  block  problem,  is 
determining  the  conditions  under  which  such  a  case  can  occur.  That  is  not  addressed 
in  this  report.  The  special  case  of  the  top  block  rocking  only  is  essentially  the  problem 
of  a  single  rocking  block  which  has  been  studied  extensively  in  the  literature  and  will 

be  considered  only  minimally  in  this  report. 

/ 

Ca^es  1  and  3  are  similar  except  for  the  angular  rotation  direction  and  the  same 
can  be  said  for  cases  2  and  4  in  terms  of  their  similarity.  What  this  implies  is  that  cases 
1  and  3  will  have  the  same  form  of  the  equations  of  motion  except  for  a  sign  difference 
on  the  RHS  of  the  equation,  just  like  the  single  rocking  block  problem.  Case  2  and  4 
will  exhibit  the  same  characteristic,  thus  lending  some  efficiency  to  representing  these 
equations  in  matrix  form.  Applying  a  Newtonian  formulation  to  these  models,  the 
specific  method  used  to  derive  the  equations  of  motion  is  rather  straightforward  and  is 
presented  next. 

To  formulate  two  independent  equations  of  motion  for  each  of  the  four  cases 
under  consideration  in  this  report  the  dynamic  moments  were  equated  to  the  sum  of 
the  moment  of  the  forces.  More  specifically,  it  is  possible  to  derive  three  equations  of 
motion  for  this  model,  but  only  two  equations  are  independent.  There  are  also  internal 
forces  in  this  system  due  to  the  block-to-block  contact  and  spring  force  that  complicate 
the  derivation.  To  simplify  the  derivation  as  much  as  possible,  the  first  equation  of 
motion  was  derived  for  the  top  block  at  its  point  of  rocking  with  respect  to  the  inertial 
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CASE  1 

e>o 

a>0 


Y 


CASES. 

0<O 

a<0 


CASE  2. 
0>O 

a<0 


Y 


CASE  4 

0<O 

a>0 


FIGURE  4  -  Four  Phases  of  Motion  Under  Consideration 


reference  frame  E.  This  alleviated  the  need  to  consider  the  blockrto-block  contact 
force  since  the  moment  of  such  a  force  about  the  point  of  application  is  always  zero. 
However,  since  the  point  of  block  B2  at  which  the  moments  are  determined  is  not  fixed 
with  respect  to  the  inertial  reference  frame,  this  must  be  accounted  for  when 
calculating  the  dynamic  moment  of  this  block  and  results  in  an  expression  that 
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includes  some  Coriolis  terms.  The  second  equation  of  motion  was  derived  for  block 
Bj  about  the  point  of  rocking  for  the  system  which,  is  coincident  with  the  point  of 
rocking  for  block  Bj.  This  alleviated  the  need  to  consider  the  system  contact  force 
since  the  moment  of  such  a  force  about  the  point  of  application  is  zero,  as  before. 
However,  the  spring  force  and  the  contact  force  due  to  the  top  block  must  be 
accounted  for  since,  in  general,  the  line  of  action  of  those  forces  will  not  pass  through 
the  point  of  application  and,  hence,  their  moments  will  not  be  zero.  This  results  in  one 
equation  of  motion  per  block  for  each  of  the  four  cases.  Take  case  1  for  example  (see 
Figure  4).  Equation  12a  is  the  result  of  taking  the  moment  of  the  forces  acting  on 
block  B i  about  point  G\  (see  Figure  3),  which  is  fixed  in  the  inertial  reference  frame  £, 
and  equating  these  forces  to  the  dynamic  moment  of  block  Bj.  The  forces  acting  on 
block  Bj  arc:  the  weight  of  block  B  j  in  gravity  field  g,  the  spring  force  B2  (always  in 
tension),  the  contact  force  at  point  Q\  due  to  block  B2,  and  the  contact  force  of  the 
system  at  point  G\.  Since  point  Q2  of  block  B2  is  coincident  with  point  (?i  of  block  Bj 
there  is  no  extension  of  spring  Bj  so  spring  Bj  does  not  contribute  to  the  dynamics  of 
this  particular  case.  Also,  since  the  moments  of  the  forces  are  being  summed  about 
point  G\  (which  is  fixed  in  the  inertial  reference  frame  e  for  this  specific  case)  the 
moment  of  the  contact  force  at  point  G\  is  zero.  Thus,  the  moments  of  the  three 
remaining  forces  about  point  G j  arc  equated  to  the  dynamic  moment  of  block  B  j  about 
point  Gj  to  derive  equation  12a.  Since  this  is  a  2DOF  system,  the  second  equation, 
equation  12b,  is  derived  by  evaluating  the  moments  of  the  forces  acting  on  block  fl2 
about  its  point  of  rocking,  point  Q2,  with  respect  to  the  inertial  reference  frame  £,  and 
equating  these  moments  with  the  dynamic  moment  of  block  fl2,  The  forces  acting  on 
block  B2  arc:  the  weight  of  block  B2  in  gravity  field  5, the  spring  force  B2,  and  the 
contact  force  at  point  Q2  due  to  block  Bj.  Since  point  Q2  of  block  B2  is  coincident 
with  point  Q\  of  block  Bj,  there  is  no  extension  of  spring  Bj,  so  spring  Bj  does  not 
contribute  to  the  dynamics  of  this  specific  case.  The  moment  of  the  contact  force  at 
point  Q2  is  zero  since  it  is  being  summed  at  the  same  point  also.  This  same  procedure 
is  applied  in  the  other  three  cases  and  the  general  equations  are  presented  next. 


CASB.i.»v>J  ZSLS^Jl 

1 Obici  +  mi(ni2  +  h2)  +  mj/ffc/-/;)2  +  (2ajf]}  9  (13a) 

-m2{[h2(brl1)-2n2a]]sin(a-Q)  -  [n2(brlj)+2h2aj]cos(a~B)y&  2 
-m2f[n2(brlj)+2h2aj]s\x[(,a~B)  +  [h2(bj-lj)-2n2cij]cos(ct—B)}(ii 
■[m1h]g+2m2a1g]swB  -  K2(b2+l2)(lj+l2)sin(a-B )  -  -[m1n1g+m2g(b14})JcosQ 

@B2C2  +  m2  ln22  +  b22]}6.  (13b) 

-m2{[h2(brl])-2n2aI]sin(a-B)  -  [n2(brlj)  +  2h2aj]cos( a-0)>02 
+nt2'([n2(bi~h)'^2h2Qi]sin(CL—B)  +  [h2(b  j-lj  )-2n2a  j]  cos(cl-Q))B 
=  • m2g[n2cosa  -  h2sina]  -  K2(lj+l2)(b2+l2)sin(Oi-B ) 
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CASE2->e>Q,g<e 


(Ibici  +  ml(nf )  +  m2[(b]+l]f+(2ai)?]}Q  (14a) 

mtn2{[h2(b]+l2)+2(2b2mn2)ailsin(Q’~Qi)+[(2b2mn2)(b]+l2)’2h2Cii]cos(ct—Q)yct  2 
+m2{[(2b2-n2)(bj+l2)-2h2a]]sin(a-Q)-[h2(bj+l2)+2(2b2-n2)ai]cos(a-Q)}Ct 
-KI(b2+l])(li+l2)sin(a-Q)-(2m2gai+m1hl8)sinQ^  -[m]b]g+m2g(b1+l2)]cose 

"Ob2C2  "l"  rn2[(2b2mn2)n2  4“  b22  ]}  (5t  (14b) 

-m2{[h2(bi+l2)+2(2b2-n2)ailsin(a-Q)+[(2b2-n2Xbi+l2)-2h2aj]cos(  a-9)}§ 
-m2{[2b2- n2Xb]  +l2)m2h2a  ]]sin(o,—6)-[h2(bi +l2)+2  (2  b2-n2Xl  ]]cos(a.~Q) }  0 
=  m2g[b2cosa,Jfh2sina]  -  Kj(b2+liXlj+l2)sin(a-B ) 

CASE  3  -fr.9,.*.Q . SLL& 


{lBlCl+ml[(2bi-n])P +hf]+m2[(bj-l2ji+(2aij2]}Q  t_  (15a) 

+m2{ih2(brh)m2(2b2-n2)oiJsin(a-Q)+[(2b2-n2Xbrl2)+2h2a}]cos(a-Q)y^: 
-m2{[(2b2-n2Xb]~i2)+2h2a]Jsin(a-Q)-[h2(bj-l2)-2(2b2-n2Xii]cos(a-B)}(L 
-K1(b2+liXh+h)sin(a-Q)mimlhl8Jt'2m2al8]sin&  =  [mjg(2b1-n1)+m2g(b}-l2)]cosQ 

{lB2C2+m2l(2b2~n2)i +h22]Y&  ..  (15b) 

+m2{[h2(b]-I2)'2(2b2-n2)ai]sin(a-Q)+[(2b2-n2Xbrl2)+2h2a]]cos(a-&)}% 

+m2{[(2b2-n2Xbrl2)+2h2a1]sin(a-Q)-[h2(brl2)-2(2b2-n2)a1]cos(a-Q)}Q 
■-  m2g[(2b2‘n2)cos<X+/i2*inCi]  -  ^/(b2+ljXIj+l2)sin(a-0) 

CASE  4  6  <;  Q . 

{iBlCl+Mltftbi-ftif+hfj+^ffbi+lif+fiajf  ]}Q  "  (16a) 

+m2{[h2(bj+l1)+2n2ai]sin(a-9)-[n2(bi+lj)‘2h2a1]cos(a~B)ycL  2 
+m2{[n2(bJ+l1)-2h2ai]sin(a-9)+[h2(b1+li)+2n2a1]cos(a~B)}CL 
-K2(b2+l2)(li+h)sin(a-fy‘[mlhl8+2m2ajg]sind~  [mibig+m2g(bi+l})]cosQ 

{lR?C2+m2(n22+h22)y&  (16b) 

+nt2{lh2(bj-¥lj)+2n2aj]5in(a.-9)-[n2(bj+lj)-2h2aj]cos(a-B)}d2 
-m2{[n2(bj+l])-2h2a}]sin(ct.--B)+[h2(b]+li)-¥2n2aj]cos(CL~B)}Q 

-  nt2g[h2sina-n2Cosa]-K2(b2+l2Xh+h)sin( 


where  /g;ci  and  /fl2C2  arc  mass  moments  of  inertia  of  each  block  about  their 
respective  centers  of  mass.  Note  that  equation  (a)  refers  to  the  equation  of  motion  for 
block  B\  and  equation  (b)  refers  to  the  equation  of  motion  for  block  B2>  Also  note  that 
a  single  equation  of  motion  for  the  whole  system  can  be  derived  and  would  just  consist 
of  the  sum  of  equations  (a)  and  (b),  which  is  to  be  expected. 


16 


These  are  the  general  form  of  the  equations  of  motion  and  they  are  extremely 
nonlinear  so  some  simplifying  assumptions  are  desired  to  make  these  equations 
more  tractable.  As  before,  assume  a  homogeneous  distribution  of  mass,  small  angular 
displacements  and  that  the  connecting  springs  are  located  symmetrically  about 
thevcrtical  axis  of  each  block.  The  equations  are  presented  next. 


CASE  1  9  >  0  .  a  >  9 


{lBlCl+ml(bi2+<*i2]+m2[(bimhfi+(2alP)y&  +  fn2[b2(b]-l})+2a2a]]'6.  (17a) 

■  ^2(b2J'hXh+h)(X‘  '  lmlai8+2m2a]8)~K2(b2+h)(h+h)]® 

=  ■[m1b1g+m2g(b1-l1)] 

m2[b2(b]-lj[)+2a2aj]Q  +  {tB2C2+m2fb22 +a22]y&  (17b) 

-/ m2^28 mK-2(b2 +h)(h+h)]Q-  "  K2(b2+h)(h+h)Q  -  'm2b2i 

CASE  2  -->  6  >  0  .  a  <  6 


{lBici+mi(bid'+Q2)+in2[(bj+l])*+(2aif]}§-m2[b2(b!+l2)-2a2aj]&  (18a) 

•[ntja]g+2m2aig-K](b2+l])(li+l2)]Q-Ki(b2+li)(li+l2)Q‘!=m[rnibi8+m2(bi+l2)] 

•m2[b2(b] +l2)-2a2Q]] 0  +  {lB2C2+m2(b22 +a22 )}&  ■  Kl(b2+hXh+h)Q  (18b) 

-[m2a2g~Kj(b2+liXli+l2)]a>  =  m2  b2& 

CASE., 3 

{lBici+mifbi2+a]2]+nt2[(bj-l])?+(2ajf]}Q+  m2{b2(b]-l2)+2a2aj}‘6,  (19a) 

■  lm]a]g+2m2ajg-Kj(b2+li)(li+l2)]Q  -  Kj(b2+hXh+h)a  ~  mibig+^&ibrh) 

m2[b2(bj~l2)+2a2ai] 9  +  {lB2C2+m2[b2+ay ]}CL  ■  Kjffy+ljXh+h)®  (19b) 

-  [nt2a2g~Kj(b2+li)(li+l2))0.  =  m2&2£ 

CASE  4  -->  Q  <  0 .  a  >  9 

{lB!Cl+ml[b2+a2]+m2[(b]+lij2+(2aij2]}Q  -  +h  )’^a2al]  &  (20a) 

•  [mia]8+2,n2ai8mK2(b2+h)0l'+'h)J  ^2(b2+hXh+h)^  ~fnibjg+m2g(bi+lj) 

■  m2[b2(bi+lj)-2a2aj]Q  +  {lB2C2+m2(b22 +022 *  ^2(b2+hXh+h)  0  (20b) 

-  lM2a2&'K2(b2+l2)(h+h)l  0  = 


To  better  visualize  the  similarities  between  the  equations,  they  are  presented  next 
in  matrix  form,  which  lends  a  better  understanding  of  their  symmetry,  First,  some 
simple  substitutions  will  be  made  for  conciseness. 
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Let 


h  -  hici+mdbi+ai] 

Iq  —  lB2C2+m2lb2*'*‘a2*] 

A  -  [b2(b]-l])+2a2a]]  =  [b2(bj-l2)+2a2aj] 

B  -  [b2(bi+l2)-2a2aj]  =  {b2(b]+lj)-2a2aj] 

since  l]  -I2  because  of  symmetric  spring  placement. 
*Sl  -  Kifh+hXli+h) 

Ksz  =  ^2(b2+h)(h+h) 


The  first  matrix  equation  corresponds  to  cases  1  and  3  where  the  equations  are 
the  equations  for  each  individual  block,  i.c.,  equations  (a)  and  (b).  As  mentioned 
previously,  the  only  difference  between  cases  1  and  3  is  the  sign  change  on  the  RHS. 
This  is  accounted  for  by  using  the  sgn  function,  as  shown  below.  The  same 
nomenclature  applies  to  the  second  matrix  equation,  which  represents  cases  2  and  4. 
The  resulting  equations  are  presented  next. 


CASE  1  — >  _fl  >  0.  a  >  6  and  CASE  3  -->  9  >  0 .  a  <  8 


lG+m2[(br  h?+(2aif]  'M" 

•  1 
6 

KS2  -  mjajg  -  2m2ajg  .  ks2 

0 

m2^  Jq 

ct 

+ 

■  Xs2  KS2  ‘  ^2a28 

a 

[  mjbjg  +  m2g(br  lj))sgn  (Q-a) 
[  »t2b2S  J  sgn  (d-a) 


CASE.2_--;> ,fl £ILxl& fi and  CASE 4 8 < Q , 


I  lQ+m2t (b]+  -  m2  B 

§1  -  fniajg  -  2m2ajg  -  KS1 

1 

0 

-  B  Iq 

&J  |  •  Ksi  *31  ■  m2a28 _ 1 

a 

8 

-  [  m  jbjg  +  m2g(b1+l2)  sgn  (Q~a) 

l  m2b28]  sgn  (0-a) 

In  this  form,  the  symmetry  between  the  cases  is  readily  apparent  and  facilitates  the 
understanding,  determination  and  encoding  of  a  solution  method.  ‘This  system  of  local 
linearized  equations  is  globally  nonlinear  due  to  the  discontinuities  at  transition. 
However,  valid  and  piecewise  analytical  solutions  can  be  found  in  closed  form 
between  transitions.  The  method  to  do  so  and  the  solutions  arc  presented  in  section  3. 
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3.  ANALYTIC  METHOD  AND  SOLUTIONS 


As  shown  in  the  previous  section,  the  nonlinear  equations  of  motion  (Eqs.  13 
through  16)  in  their  most  general  form,  with  the  addition  of  an  arbitrary  forcing 
function.#/),  can  be  represented  by 


[i(r\)J  {r\>  +  U  (r\)J  <A2)  +  IK]  <L(i\)}  =  [R]  {S(r\)}M  (21) 


where  the  matrices  /,  J,  K,  and  R  are  coefficient  matrices  that  are  defined  by  the 
physical  parameters  of  the  system  of  interest  and  the  vectors  L  and  S  are  discontinuous 
functions  of  the  original  coordinates  rij,  i.e.,  r|i  -  6  and  r\2-a.  The  forcing 
function,  f(t),  is  a  general  deterministic  or  probabilistic  one  and  can  include  horizontal 
and  vertical  components  of  some  ground  excitation,  or  lateral  excitation  due  to  impact, 
wind,  or  other  effects.  As  mentioned  previously,  these  equations  are  nonlinear  with 
specific  nonlinear  terms  arising  from  Coriolis  effects  and  geometrical  constraints,  that 
contribute  trigonometric  terms  to  the  vectors  L  and  S.  The  equations  arc 
nonautonomous  under  forced  vibration.  In  addition,  the  matrices  K  and/?  undergo 
characteristic  changes  based  on  the  sign  of  the  original  coordinates.  As  a  result,  there 
arc  N  equations  to  describe  N  degrees  of  freedom  with  2s  possible  solution  sets 
depending  on  the  system  parameters.  Even  after  a  linearization  process  by  assuming 
small  angles  of  rotation,  atypical  equations  remain.  Specifically,  the  resulting  stiffness 
can  be  negative  and,  by  virtue  of  the  discontinuities,  the  systems  of  equations  are  only 
piecewise  linear. 

This  equation,  in  its  present  form,  docs  not  have  any  known  analytic  solution  but 
can  be  solved  by  a  number  of  numerical  techniques.  Instead  of  solving  the  equations 
numerically,  the  equations  will  be  linearized  and  piecewise  analytic  solutions  will  then 
be  determined  for  the  free  vibration  problem. 

After  linearization  about  0  -  0  and  a  -  0,  and  setting  the  forcing  function  f(t) 
equal  to  zero,  the  system  of  equations  become 


im>  +  Mm  ={Q>  «  (22) 


where  the  vector  Q,  independent  of  time,  represents  a  geometrical  constraint  due  to  the 
linearization  process  and  not  due  to  the  general  forcing  function.  Now,  in  general, 
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these  equations  are  inertially  and  elastically  coupled.  Utilization  of  modal 
mathematical  techniques  will  allow  uncoupling  the  /  and  K  matrices  thereby  allowing 
a  solution  set  to  be  determined  for  this  system. 

Consider  a  linear  transformation  which  maps  the  original  coordinate  system  <ri> 
into  a  set  of  general  coordinates  {0}  as  shown  below.  Using 


(r)}  =  [//]  {q} 


(23) 


where  matrix  V  is  the  transformation  or  modal  matrix  and  the  vector  q  represents  the 
new  or  generalized  coordinates  allows  the  transformation  of  equation  22  to 


mmq}  +  [K][U]fq}~  { Q}  (24) 


Now,  pre-multiply  through  equation  24  by  the  transpose  of  matrix  U  which  gives 

tvfmmw  +  ivFikuvm = iu]t{q>  od 


Now  allowing  some  straightforward  substitutions  for  conciseness,  equation  25 
becomes 


[m]{‘q>  +  [k]{q>  =  <z}  (26) 

where  [m]  =  [U)TJI][U] 

[k]  =  lUfjKUV) 

{z}  =  (U)r{Q) 


llie  object  of  these  mathematical  manipulations  is  to  uncouple  the  original  system 
of  equations.  This  is  accomplished  by  what  is  more  commonly  known  as 
diagonalizing  the  original  coefficient  matrices  /  and  K  simultaneously.  With  a 
judicious  selection  of  the  transformation  matrix  U,  this  is  always  possible  if  the 
following  conditions  at  c  satisfied:  one,  the  matrix  I  is  positive  definite  and  invertible, 
and  two,  the  /and  AT  matrices  are  real  and  symmetric.  If  these  conditions  are  satisfied, 
then  there  exists  an  orthogonal  matrix  U  that  can  be  used  to  uncouple  or  diagonalize 
the  system.  This  is  true  for  the  system  of  interest  of  this  report  but  the  proof  is  beyond 
the  scope  of  this  paper  [17].  Also,  note  that  it  is  sometimes  possible  to  diagonalize 
systems  where  the  coefficient  matrices  are  not  real  and  symmetric;  i.e.,  there  may 


20 


exist  ^ansformation  matrix  U  that  will  diagonalize  the  system  but  there  is  not 
ncviessi:  «  one!  In  this  case,  ar.  "ould  use  the  inverse  of  the  transformation  matrix 
U  instead  uf  the  transpose. 

The  next  question  one  should  ask  is  how  is  the  transformation  matrix  U 
determined?  The  transformation  matrix  U  is  just  the  eigenvectors  of  the  system  of 
interest.  Thus,  one  would  pre»multiply  through  equation  22  by  the  inverse  of  matrix  / 
(which  is  where  the  requirement  for  matrix  /  being  invertible  comes  from)  and  form  a 
new  coefficient  matrix  for  the  vector  represented  by  a  new  matrix,  A. 


+  [Aim  »  <0} 

where  [A]  =  [I]'1  [K] 


(27) 


So  the  eigenvectors  of  the  matrix  A  form  the  transformation  or  modal  matrix  U  that 
were  used  in  equation  23  to  map  the  original  coordinates,  rj,  into  the  new  or 
generalized  coordinates,  q. 

Once  this  mathematical  technique  is  used,  the  resulting  system  of  equations  ore 
uncoupled  and  lend  themselves  to  solution  by  any  of  the  various  common  means  of 
solving  ondinaiy  differential  equations  such  as  Laplace  or  Fourier  methods.  Applying 
this  modal  method  to  equation  22  results  in  equations  of  the  form  of  equation  26  with 
the  explicit  equations  being: 


3  1 

M 

*  * 

«1 

0  m2 

1 » 

«2 

*z7 

*z2 


(28) 


One  solution  set  to  these  equations  is  for  kj  >  0  and  k2  >  0\ 

zl 

qj(t)  =  En  coswjt  +  Ej2  sinwit  Ty- 

*1 

(29a) 

z2 

q2(t)  =  E2i  cosw2t  +  E22  $inw2t  Tj- 

(29b) 

or,  if  kj  <Omdk2  cOthen:  ‘ 

zl 

qj(t)  =  Eu  coshwjt  +  Ej2  sinhwjt  Ty 

(30a) 
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z2 

q2(t)  -  E2I  rpshw2i  +  E22  sinhw2t  Tt- 

*2 


(30b) 


where  the  trigonometric  or  hyperbolic  functions  comprise  the  homogeneous  solution 
while  the  last  terms  are  the  particular  solution  to  equation  28.  Both  solution  sets  need 
to  be  considered  since  the  sign  of  the  term  in  each  equation  drives  which  equation  is 
used.  For  example,  if  kj  is  negative,  then  the  correct  solution  for  qj(t)  would  be 
equation  30a.  It  is  entirely  possible  that  for  the  same  system,  k2  is  positive  and  thus 
the  coirect  solution  for  q2(t)  would  be  equation  29b.  Therefore,  these  two  independent 
solutions  would  comprise  the  solution  set  to  equation  28,  the  solution  with  respect  to 
the  general  or  transformed  coordinates.  Also,  since  the  k[  terms  in  equation  28  are 
functions  of  the  original  K  matrix  via  the  modal  mapping  process,  and  since  the 
original  K  matrix  is  case  dependent,  i.e.,  the  sign  and  magnitude  of  the  terms  can 
vary  among  the  four  cases  or  regimes  of  motion,  the  character  of  the  solution  set  for 
equation  28  can  vary  whenever  the  cases  change. 

Now  the  it i  terms  govern  whether  the  form  of  the  solution  is  trigonometric  or 
hyperbolic,  but  these  stiffness  parameters  are  valid  only  in  the  generalized  coordinate 
system.  To  understand  physically  when  the  solutions  change  character,  the  value  of 
the  ki  terms  at  the  critical  point  (fy  -  0)  must  be  understood  within  the  context  of  the 
original  coordinate  system.  As  given  by  equation  26 


W  =  [UflKIlUJ  (31) 

or  explicitly 

*1  =  ttll^ll  +  W1 1 W21  (^12  +  ^2l)  +m21^22  (32a) 

*2  “  w12^11  +  K12w22^l24^2l)  +  u22^22  (32b) 


Assuming  that  the  two  original  springs  are  placed  symmetrically,  have  the  same 
stiffness  and  substituting  in  the  explicit  expressions  for  terms  from  page  18,  then  at 
kj  =  0  and  k2  -0  expressions  can  be  derived  to  determine  the  values  at  which  the 
original  spring  stiffness  causes  a  change  in  the  character  of  the  solution. 
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K ,  = 


A  +  B 

KS)’l 

(b2  +  li)  (/j  +  l2) 

fu2\ 

1—2  — 

L 

fM21^ 

1  + 

1  \ull)  . 

(33a) 


= 


A  +  B 

„  (u22  ^ 

2  +  — 

\u12)  . 

(b2  +  l2)  (l\  + l2) 

l-2f— " 
L  \w12 / 

+ 

fu22>? 
<“12/  . 

(33b) 


where 


A  «  m^a^g  and  B  = 


Thus,  the  values  of  the  original  stiffness  terms  at  which  the  solution  changes  character 
are  functions  of  the  eigenvectors  of  the  system  which  change  as  the  original  stiffness 
terms  change  in  value.  This  means  that  for  a  specific  stiffness,  there  is  a  corresponding 
critical  value  at  which  the  solution  changes  character.  If  the  original  stiffness  changes, 
then  the  critical  value  will  change  also!  The  specifics  of  this  as  it  relates  to  this 
problem  will  be  discussed  in  section  4. 

The  wj  and  w>2  terms  in  equations  29  and  30  are  just  the  eigenvalues  of  the 
system  represented  by  equation  28  and  are  determined  by  equations  34a  and  34b. 


w,  - 


\H 

m , 


w-  = 


\H 

m o 


(34a) 


(34b) 


Since  the  mi  terms  represent  a  physical  quantity,  the  mass,  they  are  always  positive; 
therefore,  the  sign  of  the  kt  terms  determines  whether  the  solution  set  to  equation  28, 
as  outlined  in  equations  29  or  30,  contain  trigonometric  or  hyperbolic  functions. 
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Hie  Eft  coefficients  in  equations  29  and  30  are  determined  from  the  transformed 
initial  conditions,  i.e.,  the  initial  velocities  and  displacements  of  each  block 
respectively.  Since  an  analytic  solution  of  equation  22  is  known  only  in  the 
generalized  coordinate  system  or  space,  it  is  necessary  to  map  the  original  initial 
conditions  to  this  space  also.  This  is  accomplished  by  the  use  of  equations  35a  and 
35b. 


{qm  =  Mivfimm 

(35a) 

<4m  =  (muvfnHi\m 

(35b) 

where  [m]  =  IU]T[I)[U]  from  equation  26,  matrix  U  is  the  modal  matrix,  matrix  /  is 
the  original  inertia  matrix  for  the  vector,  and  the  rj  and  f)  vectors  arc  the  initial 
coordinates  for  the  displacements  and  velocities,  respectively.  In  explicit  form, 
equations  35a  and  35b  arc: 


tfl<0) 

m i  0 

“n  *2i 

*11  *12 

’ll!  (0)' 

q2  (0) 

°  m2 

“12  “22 

/21  *22 

T|2(0) 

*1  (of 

wij  0 

“ll  “21 

*11  hi 

‘it,  (0)' 

<M°)_ 

°  m2 

“12  “22 

hi  hi 

1i2(°) 

(36a) 


(36b) 


Once  the  generalized  initial  conditions  are  determined  via  equations  36a  and  36b,  it  is 
relatively  simple  to  rearrange  equations  23  or  24  to  solve  for  the  Eti  coefficients.  The 
rearranged  equations  23or  24  now  become  at  t  *  0, 

Ell  «  91(0)  ±  jt  (37 a) 

*1 

E2i  =  <U(0)±~  ■  (37b) 

*2 

and  coefficients  EI2  and  E22  are  found  by  taking  the  derivative  of  equations  23  or  24 
and  then  rearranging  the  terms  to  get, 
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(37c) 


*1  (0) 


e22  - 


MO) 


(37d) 


Thus,  equations  29  and  30  can  now  be  evaluated  to  provide  the  solution  to  the  system 
equation  28. 

These  analytic  solutions,  qj(t)  and  q2(t),  are  valid  only  in  the  generalized 
coordinate  system  since  they  were  only  possible  to  obtain  in  that  space  and  must  be 
converted  back  to  the  original  coordinate  system  to  produce  the  proper  displacements, 
velocities,  and  more  importantly,  the  physical  understanding  of  the  dynamics  of  the 
problem.  Once  the  generalized  solutions  arc  obtained  as  a  function  of  time  they  arc 
transformed  back  to  the  original  coordinate  system  by  equations  38a  and  38b. 


0(7)  =  unqi(t)  +  u12q2(t)  (38a) 

a(t)  -  u21q](t)  +  u22q2(t)  (38b) 


where,  as  before,  the  uti  coefficients  are  merely  the  elements  of  the  modal  matrix  U, 


i.c., 


W11  "12 
*21  *22 


Computationally,  this  process  is  done  at  each  time  step  to  produce  the  correct 
motion  and  time  history  of  the  rocking  blocks.  The  complete  time  history  of  the 
rocking  blocks  then  becomes  the  aggregate  of  the  motion  as  evaluated  at  each  discrete 
time  step.  However,  one  must  account  for  the  discontinuities  that  occur  at  each 
impact,  which  is  represented  mathematically  as  the  transition  from  one  set  of 
equations  of  motion  with  their  corresponding  analytic  solutions  to  the  correct  or  new 
equations  of  motion  with  their  corresponding  solutions,  by  evaluating  the  local  time 
history  and  generating  new  initial  conditions  for  the  next  regime  of  motion.  As  with 
the  SDQF  system,  the  assumption  was  made  that  by  using  conservation  of  moment  of 
momentum,  a  relationship  could  be  derived  relating  the  new  initial  conditions  to  the 
previous  initial  conditions  to  account  for  the  discontinuity  at  impact.  This  procedure, 
as  discussed  in  section  2  of  this  report,  allows  the  derivation  of  equations  that  account 
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for  the  energy  dissipated  during  the  transition  from  one  regime  of  motion  to  the  next. 
However,  instead  of  one  equation,  as  was  the  case  for  the  SDOF  system,  this  method 
provides  two  equations,  as  is  to  be  expected  for  a  two  DOF  system.  As  before,  these 
equations  arc  a  function  of  the  geometry  of  the  system  and  in  addition,  they  are 
functions  of  the  mass  and  impact  velocities.  In  effect,  these  equations  bridge  the 
discontinuities  between  the  four  regimes  of  motion  and  allow  the  calculation  of  the 
new  angular  velocities  required  as  two  of  the  four  initial  conditions  for  the  next  regime 
of  motion. 

One  important  consideration  for  a  2DOF  system  is  how  to  apply  conservation  of 
angular  momentum  to  this  system.  As  mentioned  before,  two  independent  transition 
equations  are  required  to  solve  this  problem  in  closed  form,  but  there  are  three  possible 
applications  of  this  principle:  conservation  of  system  angular  momentum  or 
conservation  of  the  angular  momentum  of  each  block.  For  this  particular  system, 
conservation  of  angular  momentum  holds  true  for  the  system  and  the  top  block  only! 
The  angular  momentum  of  the  bottom  block  is  not  conserved.  The  proof  resides  in  the 
appendix.  The  transition  equations,  which  are  rather  complex  and  involved, 
corresponding  to  this  specific  model  of  a  two  DOF  system  are  presented  below  and  arc 
mathematically  equivalent  to  the  equations  first  derived  by  Psycharis  [13]. 


Eor  8  -  transition  from  case  1  to  case  4  or  case  3  to  case  2. 


(39a) 

+  dl 

(39b) 

_  hh  ~hh  .  _  h  h 
1  i^-hh  1  h  h°{ 

j\  “  h2C2  +  m2  (a22+h2) 

^2  ~  m2 1  h  (brh)  +  (2  a<  hi)] 
j$  “  'm2  /Az  (bj+h)  *  (2  aj  hi)] 

"  hiCI  +  ml  (^l2  "  b2)  +  nt2  [(2  a  if  -  (hjliXbj+li)] 
j$  m  hlCl  +  ml  ial2*‘b}2)  +  tri2 1(2  aj  f  +  (bi+lif] 


where 

and 
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For  0  -  transition  from  case  4  to  case  1  or  case  2  to  case  3,  which  is  physically  the 
reverse  of  the  previous  transition,  the  form  of  the  transition  equations  are  the  same,  but 


die  Ji  expressions  are  slightly  different  as  shown  below. 

j [  “  7b2C2  4  m2  (a2  +^2  )  (40a) 

;2  “  ~m2  /  ^2  (bi+h)  ‘  (2  <2 1  h 2)]  (40b) 

7*3  ■  n^2  & 2  (bj'h)  4  (2  aj  hi)]  (40c) 

7*4  "  hi  Cl  4  mi  (hi2  4  b2)  +  m2 1(2  a  if  -  (bj-liXbj+li)]  (40d) 

7*5  “  hi  Cl  4  ml  (al2  4  bf)  +  m2  1(2  ajf  +  (bj-lif]  (40e) 


For  a  -  transition  from  case  1  to  2  or  case  3  to  4 

(5y  =  (c2  +  ta/2)^ 

=  (c3  +  ta/3)  (t. 


(41a) 

(41b) 


where 


and 


Co  = 


JiJ4~hh 

hh  ~hh 


.  _  hh  ~hh 
2  hh  ~hh 


a 


i 


1  ... 


_  74/5  'hh 
3  hh~hh 


ids  ~hh 
hh  ~hh 


and  j  j  -  h2C2  4  ^2  (a2  ^bf) 

h  "  'w2  ^  ^2  (bj+h) "  (2  M 

73  "  h2C2  4  m2  hl+bf)  -  2  m2  &22 

7*4  -  -/»2  [b2  (brl2)  -  (2  a;  h2)] 

7*5  “  ^icy  4  mi  (°i2  4  ^;2^  +  m2 1(2  <* if  4  (bj+^f] 

* 

7g  *  ^2  72  ^2 2  4  2  aj<?2  4  bi  (bj-li)] 

h  m  hi  Cl  4  ml  (al2  4  bf)  +  mi  [(bj-lif  +  (2  dif  +  2  bifbj-li)] 
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For  a  -  transition  from  case  2  to  case  1  or  case  4  to  case  3  ,  which  is  physically  the 
reverse  of  the  previous  transition,  the  form  of  the  transition  equations  are  the  same,  but 
the  ji  expressions  are  slightly  different  as  shown  below  in  equations  42. 


h 

h 

h 

k 

h 

k 

h 


Ib2C2  +  m2  (a2  (42a) 

m2  [  i>2  (bj-li)  +  (2  a i  hi)]  (42b) 

IB2C2  +  m2(a2*+b22) -2  mibi2  (42c) 

[b2  (bj+12)  +  (2  aj  hi)]  (42d) 

l B1C1  +  ml  ( al 2  +  b2)  +  mi  [(2  aif  +  (bj-lif]  (42c) 

m2  [2  bi2  +  2  ajd2  -  bi  (bj+li)]  (42f) 


hici  +  ml  (°12  +  b2)  +  mi  [(bj+l2f  +  (2  a  if  -  2  bi  (bj+li)]  (42g) 


It  is  also  possible  to  derive  transition  equations  using  conservation  of  energy 
principles,  but  there  are  drawbacks  to  using  such  a  method.  Though  the  derivation  of 
the  expressions  relating  the  pre-impact  angular  velocities  to  the  post-impact  angular 
velocities  is  rather  straightforward  and  somewhat  less  tedious,  the  final  expressions  are 
more  complex  and  less  tractable.  The  final  expressions  are  nonlinear  because  of  the 
velocity  squared  terms,  the  velocity  coupled  terms,  and  the  trigonometric  functions  of 
the  displacement  angles.  Ignoring  these  difficulties  still  requires  one  to  arbitrarily 
select  a  “coefficient  of  restitution"  term  to  account  for  the  energy  loss  due  to  impact. 
If  the  selected  value  for  this  term  is  too  large  in  magnitude,  then  physically  this  is 
equivalent  to  adding  energy  to  the  system  which  then  undergoes  larger  and  larger 
oscillations  until  the  system  topples.  Using  the  conservation  of  angular  momentum 
method  prevents  this  by  determining  the  upper  theoretical  value  for  the  coefficient  of 
restitution  term  from  geometrical  considerations  as  stated  before.  Note  that  it  is  still 
possible  for  the  system  under  consideration  to  topple  under  the  appropriate  conditions; 
but  if  that  phenomenon  occurs,  then  it  is  a  function  of  the  initial  conditions  and  the 
geometry  of  the  problem,  not  due  to  some  lack  of  accountability  of  energy  dissipation 
during  transition. 

Another  significant  consideration  with  transition  involving  the  computational 
aspect  is  how  to  determine  which  state  the  system  lies  in  when  one  of  the  two 
coordinates  is  exactly  at  the  transition  point,  i.e.,  0  ■*  0  or  a  -  0.  Mathematically,  the 
equations  of  motion  are  undefined  at  these  points,  which  is  why  the  transition 
equations  are  needed;  but  it  is  still  necessary  to  be  able  to  determine  which  regime  of 
motion  the  system  will  enter  after  going  through  transition.  The  reason  for  this  is  that  it 
is  necessary  to  calculate  the  correct  post-impact  initial  conditions,  i.e.,  the  correct 
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values  of  0,  a,  6  and  &.  Also,  it  is  important  to  realize  that  though  this  report 
considers  the  four  main  regimes  of  motion,  as  mentioned  previously,  it  is  possible  for 
the  system  to  rock  as  a  SDOF  system  (both  blocks  rock  in  tandem),  or  it  is  possible 
that  only  the  top  block  will  rock.  The  reasons  for  not  discussing  these  cases  in  depth 
here  has  been  explained  previously  but  the  possibility  of  these  cases  occurring  must 
still  be  accounted  for  computationally  to  obtain  the  correct  system  response. 

After  each  time  step  calculation,  the  case  or  regime  of  motion  was  determined  by 
the  values  associated  with  the  0  and  a  parameters.  Transition  occurs  when  the 
current  case  differs  from  the  previous  case.  At  this  point,  since  the  pre-impact  and 
post-impact  cases  are  known,  the  transition  equations  are  used  to  recalculate  the 
current  time  step's  angular  velocities  for  the  two  new  initial  conditions  for  the  post¬ 
impact  case.  The  displacement  angle  of  the  block  undergoing  transition  was  set  equal 
to  its  transition  point,  l.e.,  6  =  0  for  a  bottom  block  impact  or  a  =  0  for  a  top  block 
impact.  The  other  displacement  angle  retains  its  previous  value;  thus,  there  are  now 
four  new  initial  conditions  for  the  post-impact  case.  These  new  initial  conditions  are 
mapped  to  a  set  of  generalized  initial  conditions  and  the  next  local  solution  for  this 
time  step  is  determined.  When  one  or  both  blocks  are  at  a  transition  point,  then  the 
angular  velocities  are  examined  to  determine  the  post-impact  case  or  regime.  Note 
that  when  the  post-impact  case  is  determined  by  the  values  of  the  angular  velocities,  it 
is  the  angular  velocities  as  determined  from  the  transition  equations  that  are  used. 

It  is  possible  that  the  velocities  can  change  sign  and  by  orders  of  magnitude 
during  transition.  The  change  of  sign  requires  that  for  computational  reasons  one 
must  consider  the  relative  velocities  of  the  blocks  to  determine  whether  they  will  rock 
in  tandem  or  separate.  In  other  words,  there  are  some  velocities  as  calculated  by  the 
transition  equations  that  are  physically  impossible  to  attain  in  order  to  maintain  the 
physical  integrity  of  the  blocks.  The  process  is  complex  and  as  an  aid  to 
understanding  the  decision  algorithm  is  provided  in  Figures  5  through  7.  For 
completeness,  cases  5, 6  and  7  were  added  to  the  four  main  cases,  as  shown  previously 
in  Figure  4.  Case  5  is  the  case  where  both  blocks  are  rocking  together,  or 
mathematically,  where  a  =  0  always.  Once  this  condition  occurs,  the  system  will 
remain  in  this  state  at  least  until  the  bottom  block  impacts.  Then  it  is  possible  for  the 
other  cases  to  occur.  Case  6  corresponds  to  the  situation  where  the  bottom  block  is 
motionless  and  the  top  block  is  rocking.  This  is  analogous  to  the 
SDOF  system  or  the  single  rocking  block  problem.  Case  7  is  simply  the  case  where  all 
system  motion  has  ceased  and  is  included  as  a  computational  aid. 

Thus,  it  is  possible  to  calculate  and  model  the  motion  of  a  2DOF  system 
representing  a  turreted  vehicle;  however,  the  primary  motivation  of  this  method  is  the 
calculation  and  prediction  of  the  acceleration  levels,  as  a  function  of  the  gross  motion 
of  the  vehicle,  at  any  point  within  the  structure.  The  solution  of  equations  17  through 
20  in  terms  of  0  and  a  allows  the  calculation  of  acceleration  levels  as  a  function  of 
these  two  parameters  and  the  geometry  of  the  system  under  scrutiny.  The  method  and 
equations  arc  presented  next. 
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Since  it  is  necessary  to  derive  the  acceleration  equations  to  determine  the  proper 
expression  for  the  equations  of  motion,  these  equations  are  just  presented  below.  See 
appendix  A,  a  sample  calculation,  for  more  detail. 

CASE  1  ->  6  >  0  .  a  >  Q 

g(px/e)  =  (ylQ  +  xltf)i  (43a) 

g(p2/e)  =  ((fe1-;1)B-2fl162)j-((i.1-i1)62  +  2o10)£ 

+  (x2 &  -  y2&2)  v  -  (*2&2  +  ;y2&)  u  (43b> 

CASE  2  ~>  9  >  0 .  a  <  9 

&{px/t)  =  (x{Q~y{&2)i-  (y{Q  +  xx&2)i  (44a) 

0(p2/e)  =  ((fc,  +  ;2)8-2a1d2)!/-((i,1  +  i2)62  +  2a!e)i 
+  ( (2b2-x2)'&-y2&2)y+  ( (2b2-x2)&2 +  y2'&)v  (44b) 

GASE3  r->(KQ.a<;9 

g(p{/E)~  ((2feJ-j;1)S-y162)^+ ((2fe1-^1)62+y10)i  (45a) 
»(p2/t)  =  ((i>1-/2)8-2a182)j- («.1-i2)02  +  2a18)£ 

+  {(2b2-x2)'6.-y2(t2)v+  ( (2i2- x2)<i2 +y2&)jf  (45b) 

CASE  4  9  Q .  a  ?>  9 

g(p{/z)  -  ((2 bl-xl)'Q-yx$2)l+  ((2 bl-x{)&2  +y ft) [  (46a) 

»(p2/e)  =  ( (ij  +/1)8-2a182)^+  ((fcj  +  /1)82'+2a18)i 

+  (x2a  -  y2(t2)v  -  (x2(t2  +  y2(t)  £  (46b) 
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where  points  p j  and  P2  are  the  arbitrary  points  in  bodies  Bj  and  B2  where  the 
acceleration  levels  need  to  be  determined.  Parameters  xj,  yj  and  x2,  y2  are  the 
coordinates  of  points  p;  and  p2  with  respect  to  the  inertial  reference  frame  e.  The 
vectors  £,  j  y  and  v  are  unit  vectors  with  the  following  relationships  to  the  inertial 
reference  frame,  e  ( X ,  Y ) . 


CASE 


i  -  cosSy  4-  sinSJ 

j  =  cos0y  -  sin0y 
** 

y  =  cosay  +  sinay 

y  =  cosay  -  sinay 

CASE  2 

j i  =  cosey  +  sinoy 

j  «s  cos0y-  sin0y 

y  -  cosay- sinay 

y  =  cosay  +  sinay 

CASE  3  -->  0  <  0  .  a  <  0 

l  -  cos0y- sinOy 

l  as  cos0y  -  sin0y 

y  =  cosay- sinay 

y  =  cosay  +  sinay 

CASE  4  ~>  Q  <0 .  a>  fl 

£  =  cos0y-  sin9y 

j  =  cosey-  sinoy 

M 

y  =  cos  ay  +  sinay 

y  as  cosay  -  sin  ay 

Modal  analysis  enables  the  decoupling  of  the  system  equations  and  the 
determination  of  analytic  solutions  as  functions  of  the  system  parameters  in  the  form 
of  coefficient  matrices.  This  allows  the  conduct  of  parameter  studies  by  adjusting  the 
coefficient  matrices,  and,  in  the  case  of  the  stiffness  matrix,  possibly  changing  the 
form  of  the  solution  to  understand  the  effect  of  varying  spring  stiffness.  The 
implications  and  results  are  presented  in  the  next  section. 
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FIGURE  5  -  Decision  TYee  for  0  <  0  • 
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FIGURE  6  -  Decision  Tree  for  0  =  0 


33 


6<o 


III 

CASE  1  CASE  5  CASE  2 


\ 

CASE  1 


1  jr 

CASE  5  CASE  2 


FIGURE  7  -  Decision  Tree  for  0  >  0 
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4.  RESULTS 


The  analytical  solutions  presented  in  section  3  were  incorporated  within  two 
computer  programs  (FORTRAN  77)  so  that  parametric  studies  of  the  effect  of  various 
values  of  the  inertia  and  stiffness  coefficient  matrices  could  be  determined.  The  first 
program  simulated  the  single  rocking  block  (SDOF)  system  while  the  second  program 
simulated  the  more  complex  stacked  rocking  block  (2DOF)  system.  The  SDOF  code 
was  run  on  a  SUN  3/50  workstation  at  32-bit  precision  and  the  2DOF  code  was  run  on 
a  CRAY  XMP/416  supercomputer  at  64-bit  single  precision  with  a  time  step  of  0.001 
sec.  The  codes  were  also  run  with  time  steps  of  0.0001  sec  with  no  significant 
differences  in  the  values  obtained.  All  the  results  presented  here  were  generated  using 
time  steps  of  0.001  sec. 

For  the  SDOF  system,  the  block  was  considered  to  be  homogeneous  with  a 
square  base,  bj=4ft  and  half-height  of  a^4ft.  Thus,  this  block  was  cubic  and  had 
dimensions  of  8x8x8  ft,  A  density  of  145  Ibs/ft3  was  assumed,  representing  concrete. 
The  block  was  displaced  0.40  radians  in  the  positive  direction  with  no  initial  angular 
velocity  and  allowed  to  rotate  In  free  vibration.  Six  values  of  the  Kt  term  were 
compared  for  their  effect  on  the  system  response.  The  results  are  presented  in  Figure 
8. 


The  trend  shown  in  Figure  8  is  rather  evident;  as  the  Kf  term  increases  in 
magnitude,  i.e.,  the  system  becomes  “stiffer,”  the  block  impacts  sooner  (at  t  =  0,15  sec 
for  K  -  5,000,000  Ibs/ft  versus  t  =  0,63  sec  for  K  =  0)  and  the  succeeding  amplitudes 
of  each  peak  are  greater  for  larger  values  of  K.  This  is  to  be  expected  since  this 
represents  typical  behavior  for  SDOF  systems  modeled  as  simple  harmonic  oscillators. 
Holding  the  initial  displacement  constant  while  increasing  the  stiffness  of  the  spring 
constant  is  equivalent  to  adding  energy  to  the  system  in  the  form  of  potential  energy  of 
the  spring.  Thus,  the  stiffer  systems  have  more  energy  and  it  is  expected  that  they 
would  impact  sooner  and  harder  than  the  more  flexible  systems.  This  is  also  why  the 
post-impact  response  is  more  energetic  for  the  stiffer  systems  since  the  energy 
dissipated  during  transition  percentage-wise  is  the  same  regardless  of  the  stiffness. 
Note  that  this  system  is  not  a  viscously  damped  harmonic  oscillator,  though  the 
governing  equations  of  both  are  the  same,  since  the  physical  model  of  the  single  block 
includes  energy  dissipation  within  the  code  to  account  for  the  ehange  in  the  pole  of 
rotation  during  transition.  An  undamped  harmonic  oscillator  would  vibrate 
indefinitely  with  no  energy  loss  while  this  system  eventually  damps  out  due  to  the 
instantaneous  angular  velocity  reduction  at  impact  from  the  transition  equation,  which 
is  a  function  of  the  system  geometry.  Note  that  for  this  specific  example  from 
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TIME  HISTORY  OF  SINGLE  BLOCK 


FIGURE  8  -  Free  Vibration  of  Single  Block  as  a  Function  of  K 


equation  11,  C  =  0.25  so  that  the  reduction  in  velocity  at  each  impact  is  75%,  while 
using  equation  11 ,  the  energy  loss  is  roughly  94%  !  It  is  important  to  note  that  the 
springy  which  can  be  thought  of  as  an  energy  storage  element  which  neither  creates  nor 
destroys  energy,  lias  no  effect  on  the  transition  of  the  system  in  terms  of  energy.  The 
most  significant  aspect  of  this  example  is  that  even  though  the  value  of  K  increases  by 
orders  of  magnitude,  the  system  response  varies  very  slowly  until  the  difference  in 
magnitude  is  around  an  order  of  5.  Thus,  it  appears  that  the  system  response,  though  a 
function  of  K,  is  a  relatively  weak  function. 

For  the  parameters  chosen  for  this  problem,  and  using  equations  4  and  5,  the 
values  at  which  the  solution  changes  character  was  determined  and  is  presented  in 
Table  1.  Mote  that  Kcr  =  149,408.0  from  equation  7.  Thus,  it  is  evident  that  the 
character  of  the  solution  changes  when  the  stiffness  of  the  spring  was  set  at  150,000 
lbs! ft;  but  it  is  interesting  to  note  in  Figure  8  that  the  response  of  the  block  did  not 
change  significantly.  It  is  somewhat  surprising  from  a  mathematical  viewpoint  that 
when  the  character  of  the  solution  changes,  the  system  response  is  relatively  benign. 
This  is  more  evidence  that  the  system  response  is  a  weak  function  of  the  stiffness. 
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TABLE  1  -  SDOF  Character  of  Solution  at  Various  Stiffness 


K  flhs/ffl _ A  fl/secl 


0 

•12.0750 

50,000 

-8.0341 

150,000 

0.0478 

300,000 

12.1707 

1,000,000 

68.7440 

5,000,000 

392.0198 

For  the  2DOF  system,  both  blocks  were  considered  to  be  homogeneous,  stacked 
symmetrically,  and  to  have  square  bases  such  that  bj-4ft  and  b2-2  ft,  The  height  of 
the  bottom  block  was  2ai-8ft  and  the  top  block  was  2a2~4  ft ,  so  both  blocks  were 
cubic  in  nature  with  the  dimensions  of  the  bottom  block  being  8x8x8  ft  while  the  top 
block  was  4x4x4  ft. 

The  first  parametric  study  was  to  consider  how  the  system  response  varied  with 
the  ratio  of  the  densities  of  the  blocks.  The  density  of  the  top  block  was  kept  constant 
at  p 2  - 145  Ibslfi3  as  the  density  of  the  bottom  block,  p;,  varied  in  multiples  of  the  top 
block’s  density  represented  as  the  ratio  P  =  Pi/p2-  The  bottom  block  was  given  an 
initial  displacement  of  0q  -  0.4  rads  with  no  initial  angular  velocity,  while  the  top 
block  was  Initially  displaced  at  cto  *  0.5  rads  with  no  initial  velocity.  The  value 
associated  with  the  parameter  K  was  kept  constant  at  zero.  The.  response  of  the  top 
block  is  presented  in  Figure  9  and  the  response  of  the  bottom  block  is  presented  in 
Figure  10. 

Considering  the  motion  history  of  the  top  block  first,  it  is  evident  that  as  the  mass 
of  the  bottom  block  increases  relative  to  the  top  block,  i.e.,  an  increasing  P  ratio,  that 
the  top  block’s  motion  becomes  less  muted  or  gradual  until  impact  of  the  bottom  block 
at  which  time  tine  additional  mass  associated  with  the  top  block  relative  to  the  bottom 
block,  i.e.,  a  smaller  P  ratio,  allows  for  a  greater  response  of  the  top  block  in  terms  of 
amplitude  and  duration  after  impact.  Since  conservation  of  angular  momentum  is  the 
principle  that  bridges  the  two  regimes  of  motion,  this  type  of  response  is  expected  and 
reasonable.  Also,  note  that  a  comparison  of  Figures  9  and  10  show  that  the  sudden 
change  in  the  slope  of  the  P  curves  in  Figure  9  correspond  to  impacts  by  the  bottom 
block.  This  is  explained  by  the  instantaneous  velocity  reductions  that  occur  for  the  top 
block  during  impact  of  the  bottom  block.  There  are  certain  conditions  that  cause  the 
top  block  to  increase  its  displacement  beyond  the  initial  displacement  due  to  the 
interaction  of  the  bottom  block  as  shown  in  Figure  9.  As  the  higher  P  ratio  shows,  tills 
phenomenon  occurs  when  the  mass  of  the  top  block  is  relatively  less  than  the  bottom 
block’s  mass.  It  is  even  possible  under  some  conditions  for  the  top  block  to  topple  due 
to  its  coupling  with  the  bottom  block,  whereas  this  behavior  is  not  seen  with  the  SDOF 
system.  This  will  be  discussed  after  considering  the  motion  of  the  bottom  block. 
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MOTION  OF  TOP  BLOCK 


FIGURE  9  -  Vibration  of  Top  Block  as  a  Function  of  the  Ratio  P 


Considering  the  motion  of  the  bottom  block,  as  seen  in  Figure  10,  it  Is  not 
surprising  to  see  that  as  the  P  ratio  decreases  the  bottom  block’s  motion  becomes  more 
gradual  and  the  first  impact  t’me  lengthens  out.  This  is  somewhat  intuitive  in  the  sense 
that  it  would  seem  reasonable  to  expect  that  motion  of  the  bottom  block  would 
quicken  as  the  mass  of  the  bottom  block  increases.  Also,  as  the  mass  of  the  bottom 
block  increases  relative  to  the  top  block,  i.e.,  an  increasing  P  ratio,  the  curves  in  Figure 
10  approach  an  asymptote  which  appears  to  coincide  with  the  K *0  curve  in  Figure  8. 
This  is  reasonable  in  that  the  effect  of  the  mass  of  the  top  block  on  the  system  response 
decreases  with  an  increasing  P  ratio.  Figure  10  also  shows  that  once  the  bottom  block 
impacts,  the  system  motion  damps  out  very  quickly.  This  phenomenon  is  thought  to 
be  caused  by  the  coupling  inherent  in  these  type  of  systems. 

The  trend  in  Figure  9  shows  that  increasing  the  P  ratio  causes  the  displacement 
of  the  top  block  to  Initially  increase  rather  than  decrease.  Physically,  as  the  mass  of 
the  bottom  block  increases,  the  restoring  moment  due  to  gravity  acting  on  the  bottom 
block  increases,  causing  the  bottom  block  to  rotate  down  faster.  This  displaces  the 
point  of  rocking  associated  with  the  top  block  causing  a  destabilizing  moment  on  the 
top  block;  hence,  its  displacement  angle  increases.  Figure  9  shows  that  this 
displacement  increase  due  to  the  mass  ratio  does  not  grow  indefinitely  but  approaches 
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MOTION  OF  BOTTOM  BLOCK 


HGURE  10  -  Vibration  of  Bottom  Block  as  a  Function  of  the  Ratio  P 


a  terminal  value  that  is  believed  to  be  a  function  of  the  physical  parameters  of  the 
system.  Thus,  if  a  stability  criterion  for  the  top  block  was  known,  then  a  comparison 
with  this  terminal  value  could  address  the  toppling  of  the  top  block. 

Since  the  physical  criteria  for  toppling  or  system  instability  for  MDOF  systems 
are  not  generally  well  defined,  no  general  conclusions  can  be  drawn  about  the  stability 
of  MDOF  systems  at  this  time.  It  is  possible  to  generate  mathematical  constraints  or 
criteria  for  stability  for  the  linear  regimes  of  motion  in  terms  of  the  physical 
parameters  of  the  problem,  but  this  does  not  address  the  issue  of  system  stability 
because  at  transition  the  dynamics  of  the  problem  change  instantaneously.  This 
discontinuous  change  has  the  potential  to  convert  a  mathematically  unbounded 
solution  to  a  stable  one  or  a  mathematically  bounded  solution  to  an  unstable  solution. 
One  straightforward  method  to  solve  this  dilemma  is  to  analyze  the  response  of  a 
specific  system  and  see  if  the  time  history  grows  unbounded. 

a 

A  cursory  examination  of  the  dynamic  stability  of  this  particular  model  in  terms 
of  how  varying  the  stiffness  affects  the  system  response  was  performed  and  the  results 
are  presented  in  Figure  11.  For  this  example,  two  stiffness  constants  were  chosen,  k~0 
lbs! ft  and  b- 50, 000  Ibs/ft.  The  densities  of  both  blocks  were  set  at  145  Ihslfi  and  the 


bottom  block  had  dimensions  of  &x&c8  ft3  while  the  top  block  was  cubic,  also  with 
dimensions  of  4x4x4  ft3.  The  initial  displacements  of  both  blocks  were  zero  with  the 
initial  velocities  of  both  blocks  as  variables  to  generate  the  stability  plot. 


The  most  important  consideration  in  generating  this  plot  was  what  stability 
criteria  to  use.  For  the  case  where  the  stiffness  was  zero,  i.c.,  no  spring,  the  top  block 
was  considered  to  topple  after  passing  45  degrees  in  displacement  which  is  the  static 
equilibrium  point  for  a  cubic,  homogeneous  block.  It  is  harder  to  justify  toppling  of 
the  bottom  block,  at  the  same  point  for  the  following  reason.  If  the  point  of  rotation  of 
the  top  block  resting  on  the  bottom  block  is  in  the  same  direction  of  displacement  of 
the  bottom  block  then  the  additional  weight  of  the  top  block  would  impart  an 
additional  destabilizing  moment  on  the  bottom  block  causing  toppling  to  occur  sooner 
than  45  degrees.  Conversely,  if  the  point  of  rotation  of  the  top  block  on  the  bottom 
block  is  on  the  furthest  edge  of  the  bottom  block,  in  terms  of  the  bottom  block’s  point 
of  rotation,  then  the  weight  of  the  top  block  would  provide  an  additional  stabilizing 
moment  so  that  toppling  of  the  bottom  block  would  occur  after  45  degrees.  For 
simplicity,  the  critical  stability  point  was  chosen  to  be  45  degrees  for  both  blocks  since 
the  primary  interest  was  to  account  for  the  effect  of  the  system  stiffness  on  system 
stability,  and  not  to  determine  the  exact  stability  state  of  the  system. 
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The  stability  criteria  for  the  case  where  the  system  stiffness  is  not  zero  is  more 
complex.  Tn  general,  the  stabilizing  and  destabilizing  moments  are  functions  of  both 
the  spring  stiffness  and  the  displacement  angles,  so  a  simple  number  criterion  could 
not  be  determined.  The  system  was  started  and  the  response  of  both  blocks  were 
examined  after  the  displacement  exceeded  45  degrees .  If  the  solution  continued  to 
grow  unbounded,  then  the  system  was  considered  unstable  with  those  initial 
conditions.  If  the  solution  became  bounded  within  5%  of  45  degrees  or  transitioned 
which  resulted  in  a  velocity  reduction  or  reversal,  then  the  system  was  considered 
stable  with  those  initial  conditions. 

Some  general  observations  oin  be  drawn  from  Figure  11.  For  this  2DOF  system, 
stability  is  a  function  of  the  stiffness  and,  as  the  stiffnesr  increases,  the  region  of  stable 
behavior  grows  larger.  The  region  of  stability  for  both  stiffness  curves  is  centered  on 
the  origin  of  the  coordinate  system  which  is  to  be  expected,  and  the  shapes  of  the 
curves,  though  not  symmetric,  have  a  rough  order  to  them. 

The  second  parametric  study  of  the  2DOF  system  was  to  examine  what  affect 
variation  of  the  spring  stiffness  had  on  the  system  response.  The  physical  dimensions 
of  the  blocks  remained  the  same  and  the  density  of  both  blocks  was  set  at  145  Ibs/ft 3. 
The  bottom  block  was  given  an  initial  positive  tilt  of  0O  -  0.4  radians  and  the  top 
block  one  of  Oq  -  0.5  radians.  Six  values  of  K  were  used  ranging  from  zero  to 
5,000,000  Ibs/ft  and  the  blocks  were  allowed  to  oscillate  free’y.  The  results  are 
presented  in  Figures  12  and  Figure  13,  respectively. 

Figure  12  illustrates  the  motion  of  the  top  block  as  a  function  of  six  different 
spring  stiffnesses.  The  general  trend  is  somewhat  similar  to  the  SDOF  case  as  shown 
in  Figure  8;  that  is,  as  the  stiffness  of  the  system  increases,  the  curves  shift  down  and 
1  jit,  leading  to  quicker  impact  times.  For  high  values  of  K,  the  curves  appear  to 
converge  to  some  asymptote  and  the  same  first  impact  time.  It  is  believed  that  this 
asymptote  is  strictly  a  function  of  the  physical  parameters  of  the  system  and  not  the 
stiffness  terms.  Also,  note  that  the  initial  slope  of  the  curves  tends  to  be  steeper  for  the 
stiffer  systems  and  then  level  off  to  the  slope  of  the  asymptote.  To  really  understand 
this  system  response,  Figure  13  must  be  examined  in  context  with  Figure  12. 

Figure  13  illustrates  that  there  appears  to  be  little  variation  in  the  response  of  the 
bottom  block  with  respect  to  the  different  K  values.  As  with  the  SDOF  system,  it 
appears  that  the  response  of  the  bottom  block  is  a  weak  function  of  the  K  parameter 
and  primarily  a  function  of  the  physical  properties  and  initial  conditions  of  the  system. 
A  comparison  of  Figure  13  with  the  K  -  0  curve  in  Figure  8  highlights  that  the  top 
block  has  some  effect  on  the  motion  of  the  bottom  block  and  for  these  specific 
conditions,  the  effect  is  relatively  small.  One  must  be  careful  in  drawing  general 
conclusions  concerning  this  point  because  the  relative  mass  ratios  of  the  blocks  will 
confound  the  analysis  of  the  effect  to  some  degree.  When  the  mass  of  the  top  block  is 
roughly  equal  to  or  less  than  the  muss  of  the  bottom  block,  then  the  first  impact  of  the 
bottom  biock  is  celayed  with  respect  to  the  SDOF  system.  Also,  the  addition  of  the 
top  block  with  its  inherent  coupling  effects  introduces  a  synergism  to  the  system  such 
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that  after  the  first  impact,  the  system  response  really  dies  down  due  to  the  energy 
dissipation  associated  with  the  transition  equations.  This  is  even  seen  for  2DOF 
systems  that  are  so  stiff  that  they  can  be  modeled  as  one  rigid  body  due  to  a  high 
spring  stiffness,  i.e.,  large  values  of  K. 


With  this  understanding,  it  can  be  seen  in  Figure  12  and  13  that  as  the  spring 
stiffness  increases,  the  system  becomes  more  rigid  in  nature  and  begins  to  respond 
similarly  to  SDOF  systems.  That  is  the  significance  of  the  asymptote  discussed  earlier 
in  this  section.  It  represents  the  physical  condition  where  both  blocks,  though  still 
separated  relative  to  each  other,  are  rocking  in  concert  as  if  they  were  a  single  block.  It 
is  evident  from  Figure  12  that  the  spring  acts  to  pull  the  blocks  closer  to  one  another, 
which  explains  the  steepness  of  the  initial  slopes  for  the  stiffer  systems.  The  blocks 
are  just  approaching  each  other  at  a  faster  rate  and  it  is  interesting  to  note  that  the 
asymptotic  curve  starts  at  about  0.40  radians  or  the  initial  displacement  for  the  bottom 
block  since,  in  this  case,  the  bottom  biock  is  eight  times  as  massive  as  the  top  block. 
Also,  the  stiffer  the  system  the  closer  the  impact  times  of  each  block  are  to  each  other. 
For  this  example  the  impact  time  for  both  blocks  for  the  stiffest  system  is 
approximately  0.74  to  0.73  sec. 
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TIME  HISTORY  OF  BOTTOM  BLOCK 


Tim*  (sec) 


FIGURE  13  -  Vibration  of  Bottom  Block  as  a  Function  of  Stiffness  K 


For  the  parameters  used  to  produce  Figures  12  and  13,  the  critical  spring  stiffness 
required  to  cause  a.  change  in  the  character  of  the  solution,  for  each  original  spring 
stiffness  was  determined  using  equation  32  and  is  presented  in  Tfcble  2. 


TABLE  2  -  2DOF  Character  of  Solution  at  Various  Stiffness 


-5,000,000 

91,239^79,585 

68,398 

-1,000,000 

3,841,418,541 

70,802 

-500,000 

1,001,312,068 

68,374 

-250,000 

271,482,930 

68,355 

-50,000 

18,927,269 

68,251 

0 

1,803,504 

68,014 

50,000 

3,084,635 

68,857 

250,000 

192,272,695 

68,445 

500,000 

842,891,734 

68,419 

1,000,000 

3,524,577,939 

68,407 

5,000,000 

91,239,579,585 

68,397 
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Tabic  2  shows  that  for  this  specific  problem,  irrespective  of  the  original  stiffness 
chosen,  the  critical  stiffness  Kj  will  always  be  positive  and  greater  than  the  original 
stiffness;  hence,  the  generalized  solution  qj  will  always  consist  of  trigonometric  terms 
only  and  never  change  character.  However,  it  is  also  interesting  that  the  critical 
stiffness  parameter,  K2,  appears  to  be  independent  of  the  original  stiffness  but  always 
positive.  The  generalized  solution,  q2,  will  change  character,  from  trigonometric 
terms  to  hyperbolic  terms,  for  K  exceeding  approximately  69,000  lbs! ft.  Thus  the 
character  of  the  general  solution  as  presented  in  equations  29  and  30 ,  will  consist  of 
trigonometric  functions  below  an  original  system  stiffness  of  69,000  lbs/ft  and  consist 
of  a  combination  of  trigonometric  and  hyperbolic  functions  above  that  value  for  this 
particular  system. 

The  final  topic  of  this  report  concerns  the  constants  used  in  the  transition 
equations  39  aid  41  that  are  strictly  functions  of  the  physical  parameters  of  the  system 
modeled.  Table  3  presents  some  values  for  these  constants  for  various  geometries  of 
stacked  blocks.  Column  2  represents  the  values  for  the  example  problem  used  here. 


Table  3  -  C  and  D  Values  as  a  Function  of  Block  Geometry 


ci 

di 


c 2 

d2 


4x6x6 

4x6x6 

4x4x4 

.8x8x8 

8x8x8 

8x8x8 

12x4x4 

12x4x4 

0.31 

0.28 

0.75 

0.63 

0.65 

1.77 

0.44 

0.63 

0.60 

0.97 

0.81 

0.73 

0,40 

0.03 

0.19 

0.27 

0.81 

0.76 

0.37 

0.54 

0.19 

0.24 

0.63 

0.46 

1.28 

1.02 

U1 

1.18 

-1.17 

-1.54 

-0.29 

-0.52 

0.72 

0.99 

0.89 

0.82 

0.28 

0.01 

0.11 

0.18 

1,17 

0.77 

0.29 

0.52 

-0.17 

0.23 

0.71 

0.48 

Note  that  the  top  dimension  represents  the  top  block  and  the  bottom  dimension  the 
bottom  block.  Also,  the  primed  constants  are  the  values  for  the  reversed  transitions  as 
outlined  on  pages  26  through  28. 

The  trend  for  MDOF  systems  in  general,  and  for  this  2DOF  system  specifically, 
is  the  more  slender  the  system  is,  the  smaller  the  energy  loss  during  transition. 
Though  the  trend  is  somewhat  confounded  due  to  the  coupling  between  the  blocks,  it 
is  readily  evident  in  constant  cj  which  governs  the  theta  transition.  Constant  dj,  which 
influences  the  top  block  transition,  shows  that  the  most  significant  factor  affecting 
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energy  loss  is  the  relative  block  sizes,  more  so  than  the  geometry  of  the  individual 
blocks.  Constants  c3,  d3  and  d3  govern  the  top  block  transition  and  no  specific 
conclusions  can  be  drawn  from  this  relatively  more  complex  transition. 

To  summarize,  although  2DOF  stacked  block  systems  in  general  behave  in 
nonintuitive  ways  as  Psycharis  [13J  showed,  the  addition  of  spring  connectors  to  such 
systems  affect  the  response  in  a  reasonable  and  predictable  manner. 
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5.  CONCLUSION 


A  dynamic  analysis  and  computer  simulation  were  performed  on  two  simple 
block  systems.  The  first  system  modeled  was  a  single  block  with  two  attached  springs 
representing  fasteners.  The  second  system  modeled  was  two  symmetrically  stacked 
blocks  connected  together  by  two  springs  and  resting  on  a  horizontal  base.  Both 
systems  were  analyzed  for  the  case  of  hce  vibrations,  no  slip,  and  perfectly  elastic 
impacts, 

The  equations  of  motion  were  derived  using  Newton’s  second  law  and  moment 
considerations.  This  problem  is  nonlinear  due  to  the  nonlinear  terms  in  the  equations, 
and  also  due  to  the  discontinuity  associated  with  the  transition  from  one  pole  of 
rotation  to  the  other  pole  of  rotation  during  impact  of  the  bottom  block  with  the  ground 
and  impact  of  the  top  block  with  the  bottom  block.  Thus  the  system  motion  was 
considered  in  terms  of  four  main  regimes  of  motion  where  the  motion  within  a  regime 
was  continuous  and  the  system  motion  was  piecewise  continuous.  A  linearization 
process  was  applied  to  the  four  main  regimes  of  motion  only  and  transition  equations 
were  developed  to  account  for  the  physical  discontinuity  associated  with  the  transition 
to  different  poles  of  rotation.  This  linearization  process  enabled  the  resulting 
equations  to  be  uncoupled  with  modal  analysis  techniques  and  analytic  solutions  to  be 
determined.  Using  this  modal  analysis  technique  resulted  in  independent,  second 
order  equations  with  inertial  and  stiffness  coefficient  matrices  comprised  of  geometric 
and  physical  parameters  only.  Thus  a  parametric  study  was  performed  in  terms  of  the 
mass  and  stiffness  variations. 

It  is  shown  that  the  addition  of  fasteners,  represented  as  springs  of  variable 
stiffness,  affect  system  response  in  reasonable  and  predictable  ways.  The  SDOF 
system  behaves  as  a  simple  harmonic  oscillator  between  impacts  and  the  addition  of  a 
viscous  damping  mechanism  is  used  to  account  for  the  energy  dissipation  due  to  the 
transition. 

It  is  also  shown  that  the  response  of  such  a  system  is  a  weak  function  of  the 
stiffness  parameter.  The  2DOF  system  also  responds  in  a  reasonable  fashion.  It  is 
shown  that  the  spring  stiffness  parameter  plays  no  role  during  transition  and  is  at  best, 
a  minimal  function  in  terms  of  system  response.  In  fact,  the  system  response  of  the 
SDOF  case  is  affected  more  by  spring  stiffness  variations  than  the  2DOF  system. 

The  methodology  used  in  this  report  is  a  viable  way  to  handle  the  extremely 
complex  dynamical  equations  developed  for  these  relatively  simple  cases.  In 
particular,  modal  procedures  enabled  analytic  solutions  to  be  determined  and 
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implemented  within  a  computer  code  to  generate  the  time  histories.  This  report  also 
illustrated  the  need  for  a  more  in-dcpth  look  at  transition  and  its  importance  to 
properly  account  for  energy  dissipation. 


A 
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APPENDIX  A  -  SAMPLE  DERIVATION  OF  EQUAOON  OF  MOTION 


To  further  illustrate  the  process,  a  sample  derivation  of  the  equation  of  motion 
for  case  1  will  be  shown. 


FIGURE  A.l  -  CASE  1 


The  equation  of  motion  for  the  top  block  was  derived  first  by  recognizing  that  the 
dynamic  moment  of  block  B2,  8  #2/ £)  *8  equal  t0 the  sum  of  the  moments  of 

the  forces  acting  on  block  B2,  with  respect  to  the  inertial  reference  frame  e. 


$(C2»^2^£)  "  {^2  (Al) 


Point  Q2  was  chosen  to  simplify  thc  derivation  since  thc  contact  forces  acting  on  block 
B2  act  through  point  Q2,  therefore  the  moment  of  those  forces  are  zero  at  point  Q2. 
Also  thc  bar  symbol  above  B2  in  the  RH3  of  equation  Al  means  forces  external  to  Bh 
thus  the  expression  {B*  — >  symbolizes  the  forces  external  to  block  B2  that  act 

on  block  B2.  Now  consider  thcLHS  of  equation  Al  first.  By  definition, 

§(G2»B2/£)  "  Jf[?(C2’^2/£)l|E  +  ^2c2Xm2^(C2/£)  (A2> 
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where  g  ( C2>  Bj/e)  is  the  moment  of  momentum  (angular  momentum)  of  block  B2 
about  point  c2  with  respect  to  inertial  reference  frame  e  and  g  (C2/e)  is  the 
acceleration  at  point  C2  with  respect  to  reference  frame  e.  The  moment  of  momentum 
is  defined  by  equation  A3  in  general  and,  for  this  particular  case,  results  in  equation 
A4  (see  Appendix  B  for  more  detail  on  moment  of  momentum). 


g(A,S/e)  -  msi;AGsXv(Qs/E)  +[S(QS>  S)^>  (S/t) 
tAQsX}y(S/£)  XmsCQsGs 


g(C2,B2/e)  -  j£^b2C2^  e  " 

The  acceleration  ft  (C2/e)  is  simply  determined  by  equation  A5  below. 
*<cye)  -  &h(C2/e)y  -  -4lL0cJ 


(A3) 

(A4) 


(A5) 


thus  0(C2/c)  -  -d~[-blX+2all  +  (bl  -  +  )>2y]  I  using 

dr  ~  le 

Figure  3  where  the  coordinates  x2  and  y2  have  their  origin  at  Q2  of  block  B2,  i.c„ 
equivalent  to  the  parameters  n2  and  h2,  respectively.  After  taking  the  derivative  twice, 
the  resulting  expression  for  acceleration  is  presented  in  equation  A6. 

g(C2/e)  -  <.(b1-li)Q-2al6i)l-«,bl-ll)62  +  2alH)i 
+  (*2  (*2&2+y2a)0 


So,  after  making  these  substitutions  into  equation  A2  and  performing  the  cross  product 
multiplication,  the  dynamic  moment  is  explicitly  presented  in  equation  A6  below. 

§({22>fl2/e)  "  [/BIC,  +  m2^2+>’2)]yt2  +  + 

..  .2 

+  m2(G  +  9  )  ( b j  -lj)  (^2-y2)  sin  (a -6)2 
+  m2(0  +  d  )  (^j  -/j)  (*2  +  y2) cos  (<*-9)2  (A6) 
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Now,  considering  the  RHS  of  equation  Al,  the  sum  of  the  moments  of  the  forces 
acting  on  block  B2  will  be  further  delineated  into  contact  forces  (surface  forces)  and 
body  forces  (in  this  case,  gravitational  fields  only). 


M(Q2,  {B2-»B2>c)  +M(Q2,  {B2->B2}s) 

(A7) 

Considering  the  first  term  on  the  RHS  of  equation  A7,  the  contact  forces  acting  on 
block  B2  consist  of  the  force  generated  by  the  contact  with  block  B}  at  point  Q2  and 
the  spring  force  R2  acting  on  block  B2  at  point  P2,  Since  the  moments  are  taken  about 
point  Q2,  the  forces  acting  through  this  point  contribute  no  moments!  The  moment  at 
point  Q2  due  to  the  spring  force  is  given  by  equation  A8. 

¥  (02»  )  =  £Q2Pi  X  ^P2P[  (A8> 

where  k2  is  the  spring  stiffness  constant  and  Cp2pt  *s  the  spring  displacement  such 

^at  Ep2P[  m  ^2  ~'^ClP2^  m  ^2^h+h^~~  +  thus’ 

M(Q2,  {S2-»B2}C)  -  -*2  (i  j  +  /2)  (<>2  +  i2)  sin  (a  -  0)  j  (A9) 


Now  the  second  term  on  the  RHS  of  equation  A7  represents  the  moment  of  the 
field  or  body  forces  acting  on  block  B2  which,  in  this  case,  is  gravity. 


¥(Q2,  {b2-*b2}8)  ~  Cq2c2><  m  m2g(y2s‘ma-x2cosa)z (A10) 


So  substituting  equations  A9  and  A 10  into  equation  A7  gives  the  correct  expression 
for  the  sum  of  the  moments  of  the  forces  acting  on  block  B?i. 


y(Qv{B2->B2})  -  -k2(l{  +  l2)  (b2  +  l2)  sin  (a  ~Q)z 

+  m2g  (^sina-jjCOscOz  (All) 


So  to  get  the  explicit  expression  of  equation  Al,  i.e.,  the  equation  of  motion  for  block 
B2  in  case  1,  equations  A6  and  All  arc  substituted  into  equation  Al  and  slightly 
rearranged  for  convenience.  The  equation  of  motion  is  presented  in  equation  A 12. 
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/£2C2  +  m2(*2  +  >’2>]a2  +  4m2(0  +  6  )a:z 

+  w2(0  +  62)  (by-ly)  (x2~y2)  sin (a-0)z 
+  m2  (0  +  &2)  (by-ly)  (x2  +  y2)  cos  (a  -  0)2 
~  -k2(ly  + 12)  (b2  + 12)  sin  (a  —  0) z^tn28  (}>2sina -x2coscl) z  (A12) 


Note  that  this  equation  differs  slightly  from  equation  13b  since  the  terms  on  the  LHS 
are  grouped  differently,  and  in  this  derivation,  the  parameters  jr;,  yh  x2  and  y2  were 
used  in  place  of  «/,  ht>  n2  and  h2,  respectively.  However,  the  equations  are 
mathetnatically  equivalent. 

Tb  find  the  equation  of  motion  for  block  Bj  the  same  procedure  is  used  realizing 
that  there  is  an  additional  contact  force  at  point  Qj  due  to  the  weight  of  block  B2 
resting  on  block  Bj. 

(A13) 


Solving  for  the  LHS  side  of  equation  A13  first  arid  by  definition, 


(Cp By/e)  ]  xntygiCy/z)  (A 1 4) 


Utilizing  the  general  expression  for  moment  of  momentum  results  in  equation  A 15. 


g(Cx,B{/t)  -  4,[7„.rSj] 


-  's.c,8* 


(A15) 


The  expression  for  the  acceleration  term  in  the  moment  of  momentum  equation  is. 


0(cye)  -^(Cj/OlL-^koc,] 


E  df 


(A  16) 


or  explicitly 


G(Cy/t)  -■42[Xyi  +  yyl\\  -  {Xy  0  ~  )>  !  £  )  [  ~  (>>  j  Q  +  Xy  )  £  (A17) 


dt 
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So,  after  some  mathematical  manipulations  the  LHS  of  equation  A13  is, 

$(G1(B,/e)  -  [lBlCt+mx(x\+y\)^i  (A18) 

Now  the  RH5  of  -  quation  A13  is  the  sum  of  the  moments  of  the  forces  acting  on  block 
Bj  and  can  be  further  categorized  as  surface  and  body  forces,  as  shown  below. 


M(0 ,,  {B,->8,})  {B1-»B1}C)+M(G1,  {B,->Bj}*) 

(A19) 


where 

M(GV  -  ^g,C,  x  miS  m  ™\S  (y^inO-A-jCOsG)  (A20) 

M(GV  -^B^}0)  "*0^X212  + CGlpxX&2  (A21) 

where  is  the  contact  force  on  block  Bi  due  to  block  B2  and  *s  the  force  due  to 

the  spring  interface.  The  spring  force  is  proportional  to  the  relative  displacements  of 
both  blocks  and  is  given  by  l  (^2  +  ^2 “  ^ k\  •  The  contact 

force  p |9  is  more  complicated  and  can  be  founefby  performing  a  force  analysis  on 
block  B2. 


f12-  ((fc1  +  /2)e-2a,62)i.  ((B1-/1)eJ  +  2fll0-*2(fc2  ,:.2))L 

+  (,r2a  ~  y2&2)  y  -  (x2&2  +  >:2^  +  *2  (^2  +  h)  )  K  •  m2^y  (A22) 


So  utilizing  the  expression  for  the  spring  force  and  the  contact  forces  (equation  A22), 
in  conjunction  with  equation  A20,  results  in  an  explicit  expression  for  equation  A19. 

M(GV  (B,  ->Bj>)  --m2t(61-/,)2+  (2 
-m2g  [  ( b j  -■/1)  cos©  -  2a1sinG]?  +  k2  (b2  + 12)  + 12)  sin  (a  -  Q)z 

-  w2<[ >2  ( bj  -  h)  -  2  x2  a;  ]  sin  (a  -  0)  -  [x2  ( b2  - 1,)  +2  yfal j cos  (a  -  6)  }a  z 

-  (  bj  -  lj)+2 y2 a/] sin  (a  —  0)  +[y2( br  lj)-2 x-arfeos  (a -6)  >&2  z 

+  mig  [>>isine -^cosGlz  (A23) 
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Now,  substituting  the  expressions  for  equations  A18  and  A23  into  equation  A 13,  gives 
the  explicit  expression  for  the  equation  of  motion  of  block  B j . 


Abici  +  mi  ( xi2  +  y?  )  +  m2 1  (h~h  )2  +  ( 2  aj  )2]  }  Q  z 

-  ^2  { [y2  ( i>i  -  li )  -  2  x2a1]  sin(<x-0)  -  [x2  ( bj ~ lj  )+2y2ctj]  cos(a-0)  }  a  z 

2 

■ m2{[x2(b1-l1)  +  2y2a1 ]  sin(a~0)  +  [  y2  (  hj  -  /; )  -  2  x2  aj  ]  cos(a-0)}&  g 
•  [  nt]  hjg  +  :i  m2  a;  g  ]  sin0  i  -  K2  ( b2  +  l2 )( lj  +  l2 )  sin(a-0)  z 
-  *  /  mi  nj  g  +  m2  g  ( bj  - 1} )  ]  cos0  g  (A24) 


56 


/ 


APPENDIX  B  -  MOMENT  OF  MOMENTUM 


The  assumption  of  conservation  of  moment  of  momentum  for  the  derivation  of 
the  transition  equations  is  necessary  for  the  analytic  treatment  of  this  problem  and  to 
avoid  using  an  arbitrary  “ coefficient  of  restitution”  to  account  for  the  energy 
dissipation  during  transition.  The  proof  is  presented  below. 

First,  derive  a  general  expression  for  the  moment  of  momentum  of  some 
arbitrary  rigid  body  5  about  an  arbitrary  point  A  that  can  have  a  velocity  relative  to  the 
inertial  reference  frame  £  as  shown  in  Figure  B.l  below.  This  is  defined  by  the 
expression 

g  (A,  S/E )  which,  by  definition,  is  given  by  equation  Bl. 

(Bl)  g(A,S/e)  =  f  CAPX}>(P/E)d;n(P) 

S 


Note  that  point  P  is  some  arbitrary  fixed  point  in  S,  as  is  point  gj,  which  is  sonu,  point 
of  body  S  where  the  inertial  properties,  I$(Qs>S),  are  known.  Point  Gs  is  the  center  of 
mass  of  body  S  and  the  total  mass  is  given  by  ms- 
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Now 


v(f/e)  =  sW|E-^W|e  +  |,[£esf]|E 

where  y  (Qs/E)  =  £  \zOQj and  £  [rQ^p\  ^  yy  ( 5/e)  X  rQsp 

so  v  (P/e)  *  yiG/e)  +i^(5/e)  xrQsP 

thus  g(A,5/e)  -  J^pX^Gs/e)  +  F(S/e)  xr^p]^m(P)  or 

s 

g  (A,  .S'/e)  - Jj^p  x  y  (G5/e)  dm  (P)  +Jt>ip  x  [y  (S/e)  x r^p]^  (P) 
5  5 

but  £^p  =  C^o  +  *GJ>  and  £AP  =  £Aes  +  £Q,P 


sc  9(A,S/e.)-\\cAG+tGsi\xy(Qs/t)dm{P) 

S 

*  f[t:Aet+ceIJ>]  *  H5765  xcQs^dm{P)ot 

s 

(B2)  g (A, s/e) -Lc  rf» (P)  X y (G5/e)  + JrG^w <p)  * X(Qs/R) 
S  s 

+  \tAQ  X  [y(S/e)  xrCsplfifw(P)  +Jr^pX  [y (S/e)  xr(2sp]rfm(P) 
5  s  s  5 

Now  for  a  conservative  system,  i.c.,  point  P  fixed  in  body  S,  the  displacement  vector 

Zq  p  does  not  vary  so  JcGjp^J(p)  Also  jzAQsdm  (P)  ~  ms^AGs 

S  s 

and  J£aqs  *  [*  (J/E)  *  caAdm  {p)  ~1aq*w  (S/E)  x  hsdm  (p) 
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since  £aQs  an^  i y  (5/e)  have  no  dependence  on  dm(P).  Now  by 

definition,  ms^Q  G  ~  Jfg  (def.  of  mass  center) 

and  the  mass  moment  of  inertia  is  defined  as  (Q^  S)  yy  (S/e)  in  equation  A3. 

(B3)  ls  (Qg,  S)w  (S/e)  =  J cQsp  X  W  (S/e)  x  CflfP]  dm  (/>) 

5 

so  making  these  substitutions  into  equation  B2  gives  the  general  equation  for  moment 
of  momentum,  as  shown  in  equation  B4  below. 

(B4)  g(A,S/e)  =  msrAGsXy(Qs/e)  +[s(Qs,S)}y(S/e) 

+  Zaqs  x  V  (s/£)  x  "‘sr&Gj 

Now,  with  this  general  equation,  the  2DOF  system  can  be  analyzed  but  the  question  is 
how  to  apply  it  to  the  system  of  interest?  At  what  points  and  which  blocks  can  the 
principle  of  conservation  of  moment  of  momentum  be  applied?  Realizing  that  two 
independent  equations  arc  required  to  solve  the  transition  problem  in  closed  form  and 
that  there  are  three  possibilities  -  one  equation  for  each  block  and  the  system  equation 
-  the  following  proof  is  offered  to  show  where  and  how  conservation  of  moment  of 
momentum  can  be  applied  to  this  particular  system. 

Consider  Figures  B.2  and  B.3  below,  which  represent  a  transition  from  case  1  to 
case  4,  i.e.,  a  0  -  transition. 


FIGURE  B.2  -  CASE  1  FIGURE  B.3-  CASE  4 
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Definitions:  a  ••  Moment  of  Momentum  (angular  momentum) 

£  -  Inertial  Reference  Frame 

E  -  System  consisting  of  block  B  j  and  block  B2  in  inertial  reference 
frame  £  ,  i.e.,  X  «  B  j  U  B2 
Cj  -  Center  of  mass  of  block  B  j 
C2  -  Center  of  mass  of  block  Bi 
tj  -  Time  immediately  prior  to  transition 
tf  -  Time  immediately  after  transition 


First  consider  conservation  of  angular  momentum  of  the  system  about  point  G2  during 
transition  from  case  1  to  case  4.  For  this  assumption  to  be  valid,  the  following 
condition  must  be  true. 


(B5)  g  (G2,  E/e)  -  constant  as  9->0, 


which  is  mathematically  equivalent  to  equation  B6  below. 


m 


(G2,X/£) 


£ 


=  0 


but  ^~g(G2,  Z/e)  =  8(G2,  E/e)  where  §  (G2,  E/e)  represents  the 


dynamic  moment  of  the  system  £  about  point  C2  with  respect  to  reference  frame  E. 
Now  the  dynamic  moment  §  (G2,  E/e)  is  equal  to  the  sum  of  the  moments  about 
point  G2  of  all  the  external  forces  acting  on  system  Z  as  shown  in  equation  B7. 


(B7)  §  (G2>  E/e)  =  M  (G2,  {E  ->  £}  )  thus 


(B8)  £[g (G2,  E/e)]]  =  J/ (Gj,  {£  ->  E} ) 


and  integrating  this  expression  from  to  tf  to  get  an  equation  satisfying  the  condition 
given  in  equation  B5  results  in  equation  B9, 
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(B9)  9  (G2,  L/e)  =  Jjjf  ?2,  {E  — >  £}  )  *//  -  constant 

h 

But  the  moments  of  the  external  forces  can  be  further  delineated  into  the  contact  forces 
acting  on  the  system  and  the  body  forces  (gravity)  as  shown  in  equation  BIO. 


(BIO  )H(G2,  {X-»2})  =  M(G 2,  {£->£}*)  +at(G.2,  {£->£}*) 


Now  since  all  the  contact  forces  can  be  considered  during  impact  to  act  through  point 
G2,  then  the  moment  of  those  forces  about  point  G2  are  obviously  zero.  Also  during 
impact,  i.e.,  as  tf-  ^  — >  0,  all  of  the  body  forces  (gravity)  are  not  varying,  thus 

tf 

4/(G2,  {£->£}c)  =  0  and  JV(G2>  {£-*£}  ■<  constant  so 

j&f(G2,  {E->E})<if  =  J^(G2,  {E~>E}*)  -  constant  and 

U  h 

conservation  of  moment  of  momentum  of  the  system  about  point  G2  is  satisfied! 

Now  consider  conservation  of  angular  momentum  for  block  B2  about  the  pole  of 
rotation,  point  Q2.  The  angular  momentum  of  block  B2  about  point  Q2  with  respect  to 
the  inertial  reference  frame  e  is  given  in  equation  B 1 1  and  comes  from  the  application 
of  equation  B4  to  this  specific  case. 

(B 1 1)  g(Q2,B2/t)  “  m2^Q2C2  *  ^ 

+  l(CvB2)y(B2/t)  +  CQjC, * V ( B2/e)  x m2iCicr 

Using  equation  B4  it  is  easy  to  show  that: 

(B12)  g(C2,B2/£)  =  W2^C,C,  *  V (G2/^)  +  l(C2>  #2)^ 

2  2 

*rCicXy{B2/t)  x»2rCA 
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Realizing  that  q  =0,  and  combining  equations  B1 1  and  R12,  gives 
2  2 

(B13)  g(Q2,B2/e)  =  m2cQiCiXy(C2/z)  -1 -g(C2>B2/e) 


For  angular  momentum  of  block  B2  about  point  Q2  to  be  conserved,  the  following 
condition  must  be  true. 

(B14)  0  (Q2,  B2/e)  -  constant  during  transition 

Using  the  same  derivation  as  for  the  conservation  of  angular  momentum  of  the  system, 
the  following  equation  is  found. 

(B15)  2f[g  (Q2>  B2 y  )  integrating 

(B16)  jp  (Q2iB2/z)\^  =  Ja/(G2>  {&2~*B2  ))^  and  as  before 

h 

M(Q2,  (b2^b2))  =  M«22,  iBi->B2}c)+M(Q2,  {B2-*b2}8) 

and  since  all  of  the  reaction  forces  produced  during  transition  pass  through  point  Q2, 
the  moment  of  these  forces  at  point  Q2  is  zero.  Thus 

M(Q2,  {B2-*  B2})  =  M{Q2>  {B2-*  B2}8)  only  so  equation  A1 6  is 

(r 

(B17)  p(Q2,B2/£)  =  jiHQ2,  (B2^B2)‘)dt 

<i 

but  as  jy- 1[  — >  0  for  transition  the  body  force  (gravity)  does  not  vary  so 

(B18)  g(Q2,B2/i)  =  J>(e2,  {S2->fl2}®  )  dt  -  constant 

f, 

and  conservation  of  angular  momentum  of  block  B2  about  point  Q2  is  satisfied. 
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Consider  block  S;  and  its  angular  momentum  about  point  G2,  as  was  done  for 
the  system’s  angular  momentum  derivation.  The  system  moment  of  momentum  is 
equal  to  the  sum  of  the  moments  of  momentum  of  each  of  the  blocks  that  comprise  the 
system. 


(B19)  g(G2,£/e)  =  g(G2,B^/e)  +  g(G2,  #2/e) 


From  equation  A5,  g  (G2,  L/e)  -  constant,  so  using  the  A  symbol  to  represent  the 
change  in  angular  momentum  as  tj—  tj  —)  0  allows  the  following  derivation. 

(B20)  Ag(G2,  S/e)  -  Ag(G2,  Bj/e)  +  Ag(G2,£2/e)  =0 


since  the  system  angular  momentum  is  conserved.  Using  equation  B4  at  point  G2, 
equation  B21  can  be  derived  after  some  manipulation. 

(B21)  Ag(G2,B2/e)  =  A [g{QvB2/z)  +m2tG^Qxy{C2/z)]  or 

(B22)  A g(G2,B2/t)  =  Ag(Q2,B2/e)  +  m2cG^xAy(C2/E) 

but  Act  ( Q2 ,  B2/e)  “  0  from  equation  B18;  therefore,  substituting  equation  B22 

into  equation  B20  results  in  the  correct  expression  for  the  change  in  angular 
momentum  of  block  Bj  about  point  G2  during  transition. 

(B23)  Act  (G2,  “  ““^2^“GiQ2  ^  (^2^^ 


Thus  thete  is  no  conservation  of  angular  momentum  about  point  G2  for  block  B}\ 


* 
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